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Summary—The rate of reaction in an industrial, high pressure ammonia synthesis reactor is expressed in 
terms of the kinetic equation formulated by Temkin and Pyznev. An estimate of catalyst activity is 
derived, affording a means of comparing different catalysts, and of measuring catalyst deterioration. 
Optimum temperatures giving maximum rates of reaction are estimated for certain conditions of catalyst 
activity and initial gas composition. 
A reasonable approach to optimum conditions is calculated for a particular design of catalyst bed, fitted 
with a means of removing some of the heat of reaction. Incorporation of the design in a plant reactor showed 
that the best operating conditions conformed closely to those calculated. 


Résumé 


haute pression, est exprimée sous la forme de l’équation conétique de Temkin et Pyzuev. On obtient 


La vitesse de réaction dans un convertisseur industriel pour la synthése de l’'ammoniac sous 


une estimation de l’activité catalytique, fournissant un moyen de comparer divers catalyseurs et de 


mesurer la baisse d’activité d'un catalyseur. 


On obtient une valeur approximative des températures optima, donnant les plus grandes vitesses 


de réaction, pour des conditions données d’activité catalytique et de composition initiale des gaz. 


Une approximation raisonnable des conditions optima a été caleulée pour une disposition particuliére 


du lit de catalyseur, comportant un moyen d’évacuer une partie des calories fournies par la réaction. 


L’ introduction de ce dispositif dans un convertisseur industriel a montré que les meilleures conditions 


opératoires coincident fort bien avec celles prévues par le calcul. 


INTRODUCTION 
Temkin and Pyznev [1] formulated a kinetic equation 
for ammonia synthesis assuming that the rate determ- 
ining step was the process of activated adsorption of 
nitrogen. A formula for the latter was postulated on 
a semi-empirical basis, and the equation for ammonia 
synthesis in a static system given as: 
Px, . Pi, 


Pxu, 


dPXu, aN, 
dt dt 


_y Pu, 
2 Pie 


=k, (1) 


where 

Ny, = no. of mols. of N,, per unit vol. of catalyst 
after time of contact ¢ hours. 

Pou,» Py. Py, ° 
H, after time ¢. 

k, and k, are the reaction velocity constants for 


partial pressures of NH,, Ny, 


ammonia synthesis and decomposition, respectively. 
At pressures above atmospheric, the perfect gas 
laws are assumed to hold. &, and k, are related to the 


r? a . 
= Kt. As K, is 


equilibrium constant A,, thus: 


ky 
ky 
known accurately [2] over a wide range of temperature 


* TeMKIN omitted the factor 2 between the rate of N, 
adsorption and NH, synthesis, and neglected the change, at 
constant pressure, of the total volume of gas with reaction 
[see eq. (2) and (3)). 
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and pressure, experimental reaction data can be used 
to determine values of the reaction velocity constants. 
k, is the constant usually calculated, probably because 
experimenters, such as WINTER [3], first concentrated 
on the kinetics of ammonia decomposition. k, should 
obey the Arrhenius equation : 


ky = by e~ Etec! RT 
where 
E 


composition cals/mol N, reacting. 


dec = &pparent activation energy for ammonia de- 
R = gas constant. 
T = absolute temperature °K. 
b, = frequency factor. 

Temkin and Pyzuev [1] and Emmerr and Kum- 
MER [4] have applied the kinetic equation to experi- 
mental data at pressures up to 100 atm for various 
space velocities, H, N, ratios and degrees of approach 
to equilibrium. The data were obtained at three 
temperatures, 370, 400 and 450° C. The equation was 
found to interpret the data satisfactorily, except that 
k, decreased with increasing pressure, and some varia- 
tion: 40000 to 55000-—was observed in Ey,.**. 


** See Tables 4 and 5. 
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The present paper describes the fitting of the 
equation to approximate plant data, to find a measure 


of catalyst activity and a means of estimating the best 


practical reaction conditions. The data were obtained 


at 245 and 300 atm pressure, and covered a temperature 


range of 370 to 550°C. 


EXPERIMENTAL 


The plant reactor, from which the data at 245 atm 
were obtained, was an upright cylindrical forging, 
containing a single vertical catalyst bed. The basket 
holding the catalyst wa: lagged and measurements 
showed the heat loss from the converter to be only 
the total heat The catalyst 


bed was therefore considered as adiabatic 


about 2% of evolved 


Gas at 245 atm pressure and of constant (initial) 
gas com position flowed continuously into the reactor. 
After being heated to reaction temperature the vas 
passed down over the catalyst bed, where it was 
partially converted to ammonia, with a rise in tempera- 
to the heat Conditions were 


ture due of reaction. 


steady, and the pressure constant, so that at each 
level in the catalyst bed, there was a constant tempe- 
rature and also a constant ammonia concentration, 
both quantities increasing with the amount of catalyst 
traversed 

After leaving the catalyst, the gas gave up some 
of its heat to the incoming gas, before flowing out of 


the reactor. 


The variables measured were the pressure, the inlet 
gas rate and composition , the concentration of am- 
monia in the exit gas, and temperatures taken at short 
intervals through the catalyst bed by a thermocouple 
moved in a vertical sheath near the 


“ hich could be 


axis of the reactor. 

Except for a more comprehensive exploration of 
temperature in the catalyst bed, the measurements 
were the normal plant readings, and no great accuracy 


for them is claimed. 


The data at 300 atm were obtained from a reactor 
similar to the above, but containing several adiabatic 
beds in series, between each of which was provision 
for adding cold unconverted gas. The measurements 
used were the pressure, the initial gas composition, 
the amounts and temperatures of the gas added be- 
fore each bed, the total amount of ammonia made 
in the reactor, and the temperatures at the exit of 


each bed. 


Chemical 
Engineering Science 


TREATMENT OF Data 

As the catalyst beds were adiabatic, various simpli- 
fying assumptions could be made. Radial conduction of 
heat was taken to be negligible, so that catalyst temper- 
atures over any cross-section were considered uniform. 
Gas and catalyst temperatures were assumed equivalent 
and longitudinal conduction of heat was neglected. 

In the case of the single bed reactor, for the pur- 
poses of analysis, the catalyst bed, of total depth 
13 ft, was considered as a series of small horizontal 
sections, of depth 6° or I ft, such that there were 
temperature measurements at the beginning and end 
of each section. 

Since all the heat of reaction appeared as increased 
heat content of the reaction gas, the amount of reaction 
in each of the above sections could be found from the 
rise in temperature of the gas over the section, using 
specific heat and heat of reaction data. The calculation 
was commenced at the top section, as the rate and com- 
position of the gas were known only at the inlet of the 
bed. The result for the top section provided inlet gas 
data for the second section. Repeating the calculation 
for successive sections pave a series of amounts of am- 
monia synthesised, the sum of which should agree with 
the total ammonia made in the reactor, as given by the 
gas rate and the inlet and exit ammonia analyses. The 
agreement actually obtained was usually within 6%. 

In the case of the multi-bed reactor, each bed was 
considered as one catalyst section. The amounts of 
reaction and the gas concentrations were calculated 
by the same method as above, making allowance for 
the extra gas added before each bed. The initial reae- 
tion temperatures in the beds subsequent to the first 
were calculated from the exit temperatures of the 
preceding beds and the amounts of cooling produced 
by the added cold gas. The agreement between the 
calculated amount of ammonia made in the reactor 
and the total amount measured was usually between 
0 and 7%. 

The kinetic equation was adapted to the continuous 
flow system and became: 

dA 2\k Py, Pi. 
dw ' Psu, 


Pru, 
* Pit. 
Direction of gas flow 


‘ ‘ ‘ ‘ 


w M* 
A kg mols/hr NH, 


Entry to 
catalyst 
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where A = flow of ammonia in kg mols/hr after the 
gas has traversed w, M® of catalyst. 

The partial pressures of the gases were expressed 
in terms of total pressure P atm, and z the mol fraction 
of NH, in the gas after flowing over w, M’ of catalyst, 


and the equation was rearranged explicitly in k,, viz. 


_ tt / z(1 —bz)' dz 


ky “a 48 (1 +2)(L*(1 —bz)*—2*] dw (3) 


where, after traversing w, M® of catalyst, the gas rate 
is V, M* hr, measured at | atm and 20° C, the tempera- 
ture is T° K and the composition in mol fraction is: 


z of NH,, a(l —6z) of H,, and 


L = 
(l—b Zeq)” 


1 — bz) of Ny. 


a 
3 \ 


where z 


eq 


pressure P atm, and temperature 7° K. a and b are 


is the equilibrium mol fraction of NH, at 


constants determined from the initial gas composition. 

As integration of the equation for an adiabatic 
system, where z is a function of 7’, would be extremely 
complex, the following approximation was adopted. 
The average reaction velocity over each catalyst sec- 
tion, given by 42 Aw, was assumed equal to the actual 
velocity dz/dw at the arithmetical average temperature 
of the section, for the concentrations of reactants and 
products corresponding to that temperature. 

The errors thus involved were, for the case of the 
single bed reactor, generally less than 1 or 2°C in 
temperature, and 0-001 in the mol fraction of am- 
monia, and were within the errors of measurement. 
In the calculation of k,, except for conditions near 
equilibrium, an error of 2° C in temperature will cause 
an error of about 9% in k,, whilst an error of 0-001 in 
the mol fraction of ammonia will affect the magnitude 
of k, by less than 10%. 
directions, the final inaccuracy of k, should not be 


As the errors act in opposite 
more than 10%. Errors of measurement may easily 
cause an additional error of the same order in the 
k, values. 

The reaction velocity constant has not been cal- 
culated for conditions very near equilibrium, for the 
errors produced by errors of measurement are then 
very large. 

Using the approximation to apply the kinetic equa- 
tion to the data of the multi-bed reactor may involve 
greater errors, because of the more scanty information 
on catalyst temperatures. It iz estimated that in the 
majority of the examples, the error incurred in the k, 
value should not be greater than 16%. 
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With regard to the poor accuracy, it is pointed out 
that the purpose of this work was to find out how much 
use could be made of the plant data as measured, 
without introducing more elaborate instrumentation 
to improve the accuracy and frequency of measure- 
ment. 


Table 1, Experimental data for the determination of k, 


Pressure: 245 atm. Temperature range: 404-555° C. Ejfi- 

ciency range: 0-09 to 0-96. 

mol fraction 
0-021 
0-674 
0-225 
0-048 
0-032 


1-000 


Representative initial gas composition : 





Gas rate 
kM hr 


Temp. Z 
C mol frac. 


z/Aw ky 


17-0 
13-1 
15-2 
19-7 
15-6 
20-4 
43-7 


52-8 


0-096 
0-095 
0-105 
0-093 
0-093 
0-066 
0-105 
0-133 
0-203 75-2 
0-138 99-8 
0-152 161 
0-209 209 
0-175 166 
0-230 191 
0-182 358 
0-137 417 
0-163 381 
0-281 562 
0-164 614 
0-136 983 
0-214 955 
0-176 1112 
0-194 1410 
0-134 1360 
0-146 2460 
0-125 2140 
0-174 2625 
0-125 2481 
0-125 3630 
0-060 3360 
0-076 5440 
0-052 11760 


0-0321 
0-0263 
0-0286 
0-0339 
0-0247 
0-0358 
0-0439 
0-0354 
0-0408 
0-0476 
00-0590 
0-0522 
0-0484 
0-0536 
0-0665 
0-0761 
0-0631 
0-0687 
0-0743 
0-0922 
0-0833 
0-O876 
0-0842 
0-0919 
0-1088 
0-0993 
0: 1050 
0-1054 
0° 1067 
0-1164 
0-1149 
0-1234 


0-09 
0-07 
0-08 
0-10 
0-07 
0-11 
0-13 
0-12 
0-13 
0-16 
0-21 
0-19 
0-18 
0-20 
0-28 
0-32 
0-27 
0-30 
0°35 
0-47 
0-43 
0-46 
0-46 
0-52 
0-64 
0-62 
0-67 
0-68 
0-73 
0-85 
0-86 
0-96 
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Tables 1 and 2 give the values of k, determined by the 
above methods. The values of log, k, are plotted with 


RESULTS 


reciprocal temperature in Figs. 1 and 2. 


Table 2. 


Pressure: 300 atm. Temperature range: 374-499" C. Ejfi- 


ciency range: 0-19 to 0-75. 


Experimental data for the determination of k, 


mol fraction 





























Representative initial gas composition: NH, 0-060 8 
H, 0-578 : 
Ny 0-192 
CH, 0-045 | 
A 0-125 
“7-000, 
Temp. z Effic. Gas rate 
c mol frac. z/Zeq Az/Aw kM* hr Ky 
374 0-075 0-19 0-008 46-2 1-55 VOL. 
376 0-131 0-34 0-008 54-1 4-46 i -_ —— l 
379 «= 120 0-31 0-011 510 5-68 oa an 1952 
385 0-081 0-22 O-O15S 46-0 507 
387 0-078 0-21 0-016 42-2 4-82 Fig. 1. Log, k, versus 1/T for five batches of new cata- 
301 0-145 0-40 0-007 54-5 7-84 lyst. Different symbols ii different batches of 
395 0-088 0-25 0-016 45-7 7-30 ee 
395 0-075 0-21 0-013 50-4 6-10 a 
396 0-107 0-30 0-016 45-1 9-08 
396 0-087 0-25 O-O15 51-8 S48 
397 0-074 0-21 0-021 32:2 5-62 
403 0-152 0-45 0-008 67-4 16-1 
404-074 0-22 0-017 33-4 5-97 : 
405 0-157 0-47 0-006 66-6 13-1 
410 0-082 0-25 0-049 32-0 22-0 : 
412 0-079 0-25 0-049 33-3 23-8 - 
413 0-100 0-31 0-014 51-9 15-4 
415 0-148 0-47 0-014 61-6 36-7 
: . 7 are ; 14-0 \ 
420 O-113 0-37 O-O15 52-7 26-0 x 
421 0-090 0-30 0-054 31-7 38-2 Fy 
422 0-143 0-47 O-O15 60-1 43-4 ov 
425 0-123 0-41 0-021 54-9 49-1 
430 0-134 0-46 O-O1S 55-1 47°3 \ 
432 0-097 0-34 0-054 31-6 59-0 20 
436 0-164 0-59 0-010 66-5 67-1 
446 0-161 0-62 0-010 5-3 87-6 
448 0-099 0-38 0-034 33-2 60-7 0 
451 0-109 0-43 0-082 31-3 174 
463 0-147 0-63 0-027 58-1 292 
471 0-124 0-56 0-071 30-8 332 2 
487 131 0-65 0-072 31-8 60: _ need - road saad 
496 0-142 0-74 0-083 30-5 1120 
499 0-140 0-75 0-062 31-6 893 Fig. 2. Loge &, versus J/T for one batch of catalyst 


operating at various temperatures and efficiencies. 
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All the data applied to new catalyst, which was 
multi-promoted iron oxide. Poisons (oxygen contain- 
ing compounds and sulphur) were probably present in 
very small quantities in every example. Measurements 
of them were not undertaken, but their concentration 
is regarded as being fairly steady between the experi- 
mental runs. The actual small variations which oc- 
curred will increase the so-called experimental error. 

Table 1 and Fig. 1 include data at 245 atm from 
five charges of catalyst. Initial gas composition was 
very similar from one run to another, and as there 
was not much variation in initial reaction temperature, 
the adiabatic nature of the reaction has resulted in 
some interdependence of temperature and efficiency, 
where efficiency = zz 


eq* 

Fig. 2 and Table 2 present the data at 300 atm press- 
ure from a single batch of catalyst spread between 
the series of adiabatic beds. Although initial gas 
composition was again fairly constant, it was possible 
to obtain some variation of temperature which was 
independent from the variation in efficiency, thus pro- 
viding a more comprehensive test of the Temkin 
equation. 

In Table 3 are summarised values of k, and £,,. 


obtained from the two sets of data by calculating the 


best straight line for the relation log, k, versus 1/7’. 


Table 3 





Press ky at 
Data 
atm 


420°C 470°C 


Table | 
Table 2. 


33-9 318 
30-1 312 


45800 
47900 





Discussion OF RESULTS 
Fig. 2 and Table 2 show that, over the experimental 
range of ammonia concentrations, efficiencies and total 
rates of flow, the Temkin equation gave values of the 
reaction velocity constant, k,, which conformed well 
with the Arrhenius equation. 

Fig. 1 and Table 1 prove that reproducible results 
could be obtained with different batches of catalyst, 
and show that the Arrhenius equation was obeyed 
under different sets of reaction conditions and wider 
ranges of efficiencies than in the case of Table 2. 

The 
300 atm are similar. 
the k, values is approximately that expected, if the 
empirical relation between the reaction velocity con- 


values of EZ, obtained at 245 atm and 
The small difference between 
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stant and pressure, mentioned in the subsequent sec- 
tion, is correct, and the concentration of poisons in the 
reaction gas is similar for the two sets of measurements. 

The experimental scatter in the data is consider- 
able, amounting to a coefficient of variation in k, of 
29° for the first set of data, and 31% for the second. 

Nevertheless, the results are good enough to de- 
monstrate that the Temkin equation in conjunction 
with the Arrhenius equation give a good interpretation 
of the ammonia synthesis reaction in the plant reac- 
tors, making possible the measurement of catalyst 
activity in terms of the reaction velocity constant k, , 
and the apparent activation energy for ammonia de- 


, 


composition: Ey... 

The large experimental scatter will make impos- 
sible the detection of real small differences in the reac- 
tion velocity constant, and will render approximate any 
calculations on reactor operation. 


AGREEMENT OF k, AND Ey... WITH OTHER PUBLISHED 
ReEsuLtTs 
The values for E 


for ammonia decomposition, accord well with estim- 


dec: the apparent activation energy 
ates found by other workers (see Table 4). 

As the reaction velocity constant varies with press- 
ure, it is difficult to compare our results at 245 and 
300 atm with others at pressures of 100 atm and less, 
(see Table 5). Moreover, as k, is very sensitive to the 
concentration of poisons in the reaction gas, it is not 
justifiable to compare different estimates without 
a knowledge of the gas purity. 

Nevertheless it is interesting to note that our results 
can be correlated fairly well with those of Emmett and 
Table 4. 


Comparison of values for Ey... from various sources 





Pressu re 


atm E dec 


Source of data References 


TeMKIN and PyzHEv 40000 {1} 
TEMKIN and PyzHeEvs calcu- 

lations on WINTERS re- 
46500 [1] and [3} 
43500 = [1] and [5] 
46600 [1] and [5] 
45.000 [4] 
46.500 [4] 
48 900 (4) 
53000 [4] 
45800 
47 900 


Ges 6 *s at me we l 
Larson and Tour... . 10 
31-6 
EMMETT 33-3 
EMMETT 66-6 
0 ea ee ee 100 
EMMETT 100 
LC. 245 
IC. 300 


LARSON and Tour . 
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Table 5. Comparison of k, 
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values from various sources 





Pressure (atm) 





some Source of data 
-_ 300 245 100 66-6 33-3 31-6 10 
450 ies — 127 135 
Larson and Tour . 424-513 781-915 714-915 
EMMETT . 196-212 237-292 346-417 
420 JI.C.1.. ; 30-1 33-9 
Larson and Tour 223-268 246 
400 TCU. . 10-7 12-5 
EMMETT . 15-2-21-6 24-1-30-8 37-2-47-4 
370 | I.C.1. . 2-02 
EMMETT . 1-27-2-91 3-19-4-95 7-33-8-95 
Larson and Tour, by an empirical relationship taining a heavier concentration of poisons than that 


k,xP-*®, which is very similar to that found by 
oMMETT: k, «x P~®® (see Fig. 3). 




















j 40 So 60 


10Q¢ P—e 


Fig. 3. Variation of k, with pressure. 


Use or k, aS A Measure or Catatyst Activiry 


Two examples are mentioned of this use of the reaction 
velocity constant. Fig. 4 presents data for five charges 
of new catalyst which operated in a gas system con- 


1 


appertaining to the data for Fig. 1. Although similar 


catalyst of similar age was used in both cases, the 
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Fig. 4. Effect of poisons on activity of new catalyst. 


Different symbols denote different batches of catalyst. 


activity of the catalyst subject to the greater amount 
of poisons is only about one third of that operating 
under the more favourable conditions. This rapid 
poisoning effect, at least partly due to oxygen-contain- 
ing compounds, is to be expected from the experimen- 
tal results of ALMquist and Back [6], on the poison- 
ing of ammonia catalyst by carbon monoxide and 


water-vapour. 
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The change in activation energy due to the increas- 
ed poisoning is obscured by the fact that the slope of 
the line correlating log, k, with 1/7 is influenced by 
the greater poisoning of those layers of the catalyst bed 
which operate at the lower temperatures. 

Fig. 5 demonstrates the difference in activity which 
was observed for two catalysts of different manufact- 


ure operating under similar conditions. 





¢ 











4 100 1200 
roe/y 
- _— 


Fig. 5. Comparison of two catalysts. 


CALCULATION OF OptiMuM CONDITIONS 
The most useful application of the Temkin equation is 
to calculate optimum reaction conditions, and then 
to find how nearly these may be achieved in a plant 
reactor. 

The determination of optimum reaction conditions 
is a wide problem, but in the example given below, the 
variables of catalyst activity, total pressure and initial 
gas composition were fixed. 

Maximum reaction velocity for a given gas compo- 
sition, is attained at a temperature where the differen- 
tial coefficient of reaction velocity with respect to 
temperature is zero. 

Differentiating eq. (2), under optimum conditions: 

@ jdA 

eT | dw 
_.9 | Esyn , Px, Pui 
 —— om 
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where £,.,, = apparent activation energy for am- 
monia synthesis in cal/mol N, reacting. 


The above expression gives 


ky = Eaec Pru, 
k,  Esyn” Phy, Py, 

which reduces to 

| Baec }3 


L 
Esy n ( l 


(4) 
This equation gives L, for any z. 
L, being related to the equilibrium constant, is an 


accurately known function of pressure, temperature 


$00-— 


mum fermperatu re 








0-10 


al fr . af” Wh 
Mol fraction of v1.2 


Fig. 6. Relation between optimum temperature 
and ammonia concentration. 


and gas composition. Hence the optimum tempera- 
For the 
example quoted, a pressure of 245 atm was chosen, 


ture 7’, corresponding to z, can be deduced. 


with an initial gas composition: 


mol fraction 
. 0-016 
. 0-686 
. 0-228 
. 0-036 
. 0-034 
1-000 





E sec 
that its value at 470° C was 325. 
was taken for £,,,,, 
E...—E£Z 


syn dec 


was taken as 47400, and the function k, such 
A value of 20800 
as determined from the relation 
= Q, the heat of reaction per mol N, 
reacting. With this information, the optimum tempe- 
rature corresponding to any mol fraction of NH, could 
be calculated. Fig. 6 gives the relation obtained be- 
tween these two quantities, and Fig.7 the corresponding 
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maximum reaction velocities. It is of interest that 
the reaction velocities, whilst they are high and de- 
crease rapidly at low ammonia concentrations, are 

low and relatively constant 





for high ammonia concen- 
trations. 

Integration of the syn- 
thesis equation makes pos- 
the 
the ideal temperature and 


sible construction of 
ammonia concentration gra- 
dients through a catalyst 
bed. The equation is inte- 
grated in a stepwise manner 
commencing with the initial 
reaction conditions of gas 
rate, composition and pres- 
sure. For a small increment 
the 
of ammonia, the optimum 





Az in mol fraction 


temperature is calculated, 
for the 
age z in the interval, from 
eq. (4). The 


arithmetical aver- 


increment of 

















0-05 0” O15 
Mol fraction of Ni, —e 


¥ 
002 


Fig. 7. Maximum reaction velocities at optimum 


conditions, 


catalyst volume 4A w is then found by substituting the 
average values of V, z, and 7’ in eq. (3) viz.: 


PrsqisyV z(1—bz)' 


48 (1+ 2)[L4(1—bz)* aH 4? 


dw 


Az is originally chosen small enough for the approxima- 
dz 


dw 
Subsequent increments 4w are calculated in a similar 


at the average conditions. 


. ‘ Az 
tion to be applied 
Iw 


manner until the whole catalyst bed has been taken 
into account. 

The resulting ideal gradients of ammonia concen- 
tration and temperature are given in Fig. 8. The unit 
of the abcissae is: volume of catalyst traversed initial 
reaction gas rate, chosen so that results can be easily 


computed for various gas rates. 


THE BEST OPERATING CONDITIONS IN A PLANT REACTOR 


Attainment of ideal conditions in a plant reactor would 
be a complex practical problem involving removal of 
the heat of reaction from the catalyst bed at a rate 
decreasing with catalyst age, and varying through the 
bed from a high figure at the beginning of the reaction 
to a low figure at the end of the reaction. Moreover, 
the first part of the bed would be operated at such a 
high temperature that the activity of the catalyst 
there would probably be quickly impaired. 

A reactor had been designed, however, to give a par- 
tial approach to the ideal conditions, by transferring 
some of the heat of reaction to partially heated gas 
which flowed in a counter-current direction to the 
reaction gas, through vertical tubes inserted in the 
catalyst bed. The rate of removal of the heat of reac- 

tion could be varied by altering 








the rate of flow of the cooling gas, 
or its initial temperature, but the 
cooling of particular sections of 
catalyst could not be varied in- 
dependently from one another. 


The 
described above were applied to 


kinetic data and methods 


forecast the output and the best 
operating conditions of this re- 
actor, by calculating the optimum 
initial synthesis temperature and 
the optimum degree of cooling. 


The same catalyst activity and 








Vol catalyst traversed / Initial gas rate —e 


Fig. 8. Ideal temperature and NH, concentration gradients. 


initial gas composition were as- 
sumed as in the calculation of the 
ideal gradients. 
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The catalyst bed was considered in small sections, 
and the kinetic equation was used to give, for each 
section, the maximum possible reaction which was 
between the heat 


consistent with a heat balance 


D. ANNABLE: Application of the Temkin kinetic equation to ammonia synthesis in large-scale reactors 


the commencement of the reaction, and the cooling 
of one section of catalyst cannot be adjusted independ- 
ently from another. These limitations still, however, 
permit of a fair correspondence between the practical 

and ideal conditions over the last 80% 





of the catalyst bed (see Fig. 9). The 
optimum 
temperatures are higher than 550° C so 
that operation at the lower temper- 


divergence occurs where the 





atures has the compensation of pre- 
serving the life of the catalyst. This 








feature, together with Fig.9 demonstrate 








Practical optimum 
OCCA! Of um 





that in future developments of the de- 
sign of reactor, the aim should be to 
improve still further the correspondence 
between practical and ideal conditions 
in the region of the peak temperature 
and subsequently. 








f.97 
C Y 
Mo. / raction of NH, ——@ 


Fig. 9. 


of reaction, the change in heat content of the re- 
action gas, and the heat transferred to the cooling 
medium. 

The best operating conditions, calculated by the 
above method, gave an output of ammonia from the 
reactor, equivalent to an exit mol fraction of 0-21, or 
a conversion: 

(exit mol fraction—inlet mol fraction) LOO, of 
19-4% (see Fig. 9). 

This compared with a maximum ammonia con- 
version of 22-0%, calculated for the ideal operating 
conditions, and a figure of 19-0% which was the maxi- 
mum conversion actually achieved in the running of 
the plant reactor. 

The agreement between the calculated and achiev- 
ed practical conditions is good, considering the large 
coefficient of variation (30% ) in the original k, values, 
from which the k, function was determined. The actual 
catalyst temperature gradient was quite close to that 
calculated, as was the amount of heat transfer from 
the bed. 

The best practical conditions also represent quite 
a reasonable approach to the ideal. The maximum 


further improvement in the design and operation of 
the catalyst bed would only increase the output by 
less than 14%, although with the present arrangement, 
the rate of heat removal is least instead of greatest at 


Comparison of practical and ideal conditions. 
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NOTATION 
= mol fraction of H, in the reaction gas, 
corresponding to zero NH, content 
kg mols/hr NH, 
constant, such that a (l1—z) is the mol 
fraction of H, in the reaction gas corres- 


of NH, 


frequency factor in the Arrhenius’ equa- 


ponding to a mol fraction z 


tion for k,, the reaction velocity constant 
for ammonia decomposition 
apparent activation energy for ammonia 
decomposition 
apparent activation energy for ammonia 
synthesis 
reaction velocity constant for ammonia 
synthesis 
- reaction velocity constant for ammonia 
decomposition 
equilibrium constant for ammonia synthesis 
Zeq 
(1 —bzeq)? 
= number of kilo mols of ammonia 
= total pressure: atm 
‘a, = partial pressures in atm of NH,, N,, and 
H, respectively 
@ = heat of reaction, cals/gm mol N, reacting 
R = gas constant cals/deg., gm mol 
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Hydroextraction IV: Radial distribution of permeability in cakes 
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Summary—Tests have been made on a 9% in. diameter hydroextractor running at 1050 and 1400 r.p.m. 
using maize starch, precipitated chalk and kieselguhr. 


Résumé— L étude d'un hydroextracteur tournant A 1050 t/m ou A 1400 t/m a été faite en utilisant de l’amidon 
de mais, du carbonate de chaux précipité et du kieselguhr. 


It has been shown [1], [2], [3], [4], [5], [6] that and the small whirling cell (Fig.2) showed hydro- 
observed flow rates, g, of wash water through hydro- extraction permeabilities within +15% of filtration 
extractor cakes follow closely the relationship permeabilities, but there remained variation within 


this 15% which appeared systematic [4]. The varia- 


; | tions of hydroextraction permeability with cake thick- 


ness, or cake weight, was not of the form expected 


























Fig. 1. The centrifuge basket with the probe arrangement. 


for streamline flow with a liquid layer of thickness as a result of the varying compressive stresses in the 
(r. — ry) within a cake of thickness (r,—r,),asshown centrifugal field. There may be a significant resist- 
in Fig. 1. The results of tests on maize starch and ance associated with the interface between packed 
chalk in both the 9in. dia. hydroextractor (Fig. 1) material and the cloth retaining the cake. 
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APPARATUS 


The 9 in. dia. hydroextractor was fitted with a probe 
mechanism of the 
either r,, or r,, at any vertical or radial position 
within The 
shown in Fig.1 incorporates a probe J on the end of 


which allowed measurement of 


the experimental range. mechanism 
rod C, which was moved vertically by the rotation of 
wheel D. 


retention of the probe-holder H in this slot allowing 


The head E carried a slotted tube G, the 


O 


o, 


Yein diam 
qo 9 

O) .o) ain 
O © 
Pe 


Yforatons 
in base 


DOA 











Fig. 2. The small cell: (a) in position for spinning, (b) in 


position for filtration tests. 


vertical movement but no rotation of (. Horizontal 
movement of the probe was caused by rotating wheel 
M, which moved the screwed tube N and rod P through 
the head L. The rotation of NV was prevented by 
a slotted guide. From scales on the moving compo- 
nents it was possible to measure the position of the 
probe point to a reproducibility of 0-Olmm. This 
accuracy is particularly desirable in the measurement 
of radii as a variation of 0-1 mm therein is important, 
perhaps affecting the estimation of permeability 
K. by 1%. 


The small cell in Fig. 2 is the complement of the 
hydroextractor in that cakes may be formed therein 
by various methods and tested in both filtration and 
hydroextraction without disturbance [4]. The tests 
0, the 
liquid level being maintained at level @ in the trans- 
parent feed tube lying on the vertical axis of rotation E 
| Fig. 2(a)]. The arrangement of the cell for filtration 


on this cell were always carried out with r;, 


permeability measurements is shown in Fig. 2(b). The 
cell was carried in the 9 in. dia. hydroextractor by an 


appropriately rigid and balanced fitting which is not 
shown in the diagram. 
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CaKE MATERIALS 


Various samples of chalk, maize starch and kieselguhr 
were tested, the results being exemplified by the data 
in Tables 1,2, and showing the differences encount- 
ered. The starch was of the type already used in 
earlier work in this laboratory. The chalk samples 
were all precipitated materials, nominally “pure,” 
The 


samples were proprietary “‘filter-aid’’ materials. 


from various commercial sources. kieselguhr 


Table 1. Thickness of cake in 9 in. dia. hydroe xtractor 


(a) Starch (alone): ry = 10-8 em 





Ww 


i] 


200 
300 
400 
500 
600 
TOO 
800 
900 
1000 
1100 


10-42 
10-24 
10-06 
9-87 
9-66 
9-45 


9-03 
8-80 
8-60 





Starch (under kieselguhr and starch layers: r, = 10-8 cm) 





. 
ro , Inner layers* 


d 


100 K 
200 K 
300 K 
300 K 


13-5 350 8 
13-5 
13-5 


27:1 


350 8S 


350.8 





(c) Starch (on kieselguhr) 





i. oe 
Ts Tv. 


cn? 
13-6 
13-6 
13-6 
8-82 19-4 
8-82 27°3 


9-83 
9-26 
9-56 
9-87 
10-26 


350 
350 
350 
500 
700 





* The “inner layer” is that within the radius of, and thus 
exerting compression upon, the layer of material under con- 
sideration. In (g) above, the cakes were formed from three 
layers, the kieselguhr being in the middle. The mass of an 
inner layer is specified as “100K” for 100g kieselguhr, 
or “3008” for 300 ¢ starch. 
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(d) Kieselguhr (alone): r, = 10-74 em EXPERIMENTAL TECHNIQUE 
" r y—r The methods used in this laboratory for the measure- 
c 0 c . 
g em cme ment of hydroextraction permeability have already 
been described in detail and assessed as comparative 
200) 10-13 13-0 techniques [3]. In the present tests, the permeabilities 
q =o 4. ° . . . 
= “— nals in the 9 in. dia. hydroextractor were obtained from 
om O-5 25-2 , , 
. “ a measurement of the time, ¢, required for the water 
aw 9-18 315 any , 
af layer within the cake radius r, to move from r, tor, . 
600) 8-82 37-0 ti Ty 1 “2 
where at all times r, <r, < 1%. ie probe me- 
700 S47 43-9 5 L, L. I 


chanism is a very sensitive detector of the liquid surf- 





(ec) Kieselguhr (under starch): ry = 10-8 em ace radius r, and the above is the preferred technique, 





measuring the movement of r, when the liquid layer 














. % ’ ag a is draining through the cake without any feed to the 
: = = . basket. The equation appropriate to this case is 
200 10-26 11-6 TOO S o.@.tom = +e 
$s |e 1. Sag — (2) 
(f) Kieselguhr (on starch): 4a* n® Ke | vo—TL, 
" “ ae Series of cakes were tested at various speeds of revo- 
i" pi per — lution of the basket, and data for representative ex- 
periments are quoted in Tables 3, 4. In all the hydro- 
100 0-45 9-1 6-4 extractor tests the cakes were formed on two layers 
200) 45 S74 12-8 of twill cloth supported by one layer of duck weave. 
300 9-45 8-38 19-1 The dimensions of the cake were found by using 
the probe to trace the form of the cake and cloth faces. 
(g) Kieselguhr (in between starch layers The cakes in the hydroextractor were produced with 
Ww re r Ay Inner layer faces vertical to within | mm by disturbing and lightly 
PF cm cm cm? q re-moulding the face during the formation period. 
This treatment has no effect on the cake permeability 
loo LOLS 9-83 tl 30S when the operation is carried ot with a submerged 
200 LoS 9-56 13 350 5S cake and suitably delicate touch. Cakes formed in this 
w0 10-15 926 170 5 way are also reproducible for any specific quantity of 
200 10-42 0-87 11-4 5008 


material. 





- ; The radii of the cakes were measured “in situ” by 
Table 2. Ejjective mean packing densities (g dry solid cc. cake) 5 





the probe mechanism as well as by measurements of 


Material Sampte | Fucking conetly samples cut from the cake. This gave a comparison 
In Gin. dia. hydroextractor basket: with the earlier technique [1], [2] wherein cakes were 
Starch . rN dee ! OTS formed in layers by successive additions of material, 
Starch . . . ake Neo 2 74 the layer interfaces being dyed and the layer thick- 
Ma a 6-6 @ 4 er l O57 nesses obtained from sections of the final cake. 
Riner fo ly). Ses ? 0-64 The cakes in the small cell were produced using 
Chalk. 2... 1 eee eee 3 0-67 a “star gauze” during the formation period and were 
Rleasiguhe alone. 2... s l wea tested with a loosely packed glass wool pad on the 
sna of comyle 5 enter j belted inner face of the cake. These precautions did not 
Ce s 4 $0.6 6 we OR . O51 a 
affect the measured permeabilities of the cakes and 
In small cell: were essential in forming and maintaining a uniform 
Starch I 0-75 cake in the cell by disrupting the helical motion of the 
Chalk . l 0-57 fluid adjacent to the inner face of the packed mat- 
Kieselguhr alone . 2 0-37 satio erial [4). 
“a of sample 2 under | hele] In the present tests the cakes were formed on four 
ee ee ee 0-41 


layers of duck weave cloth. 
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Table 3. 
10-8 cm; "Lh, 


7-64 cm; rr, 
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Starch in 9 in. dia. hydroextractor 


83lem; W, = 160¢ 





r 
ke "Le 


10-24 O-O515 


300 
400 
5OO 


0-07225 
9-85 
9-65 
45 
0-24 
9-03 


0-0925 
0-1126 
0-1338 
0-1559 
O-1791 
(2035 


600 
700 
SO) 
900) 
1000 
1100 


S-S15 


8-60 02281 110 


ti 100" 


4 ** 
K 1050 


, 
gq sec* 


K(W+ W,)** 


» 2 
q*, aec* 


** 
mean 


2-39 2-36 
2-745 2:7 

2-99 2-99 
3-17 3-14 
3-31 3°27 
3-49 3:44 
3-62 3-59 
3-74 3:77 
3°85 3-87 


1090 
1535 
1980 
2400 
2830 
3300 
3820 
4350 
4860 





* Subsripts 1050, 1400 refer to n’ r.p.m. used in the test. 


Cakes formed in the small cell were tested for 
filtration permeability, after the centrifuging tests, by 
using the arrangement of Fig. 2(b) and measuring the 
time of fall, t, of the liquid level from H, to H, in the 
The flow 


rates used were such that the flow was laminar in all 


vertical tube of cross sectional area “‘a.” 


tests. The equation for the analysis of the filtration 
data is H, Kp. At 


H, : a.gu.b (3) 


log, 
where the filter bed is of depth L and cross sectional 
area A. 


Table 4. 
10-74 em; ry, 


Kieselguhr in 9 in, dia, hydroextractor 


7-64 cm; rL, 8:3l cm; W, 


9 — 40g 





ro K** 


log, K( W-+ W,)** 
r 


q sec? g? sec? 

768 
1230 
1680 
2150 
2650 
3260 


0-0570 vO 
10-2 
IL 
13-0 
14-6 


15-9 


10-14 
9-83 
9-51 
9-18 
8-835 
8-465 


200 
0-O888 
01222 
0-1571 
01959 
0-2385 


300 
400 
SOO 
600 
700 





* n’ 1050 r.p.m, used in all tests. 

** Values in these columns have been multiplied by 10°. 

Cakes were made in both forms of apparatus by 
building cakes from layers of different materials as 
well as homogeneous cakes, in order to cover the pos- 
sibility of variation of compressible effects for each 


material in all cake positions. 


EXPERIMENTAL RESULTS 
Table 1, contains the dimensions of the cakes used in 
the hydroextractor, the measurements including starch 
and kieselguhr layers in various positions in a lam- 


** Values in these columns have been multiplied by 10’. 


inated cake. The data show the variations of radii 
in the apparatus and the reproducibility obtained. 
The data are plotted in Fig. 3 and show that to the 


accuracy of the measurements the packing density of 














6 04 


? 
< 
—@& 


Mass x 10 


Fig. 3. Relation between volume of cake in hydroextrac- 
tor and dry weight of material in it: A Starch alone, or 
starch on or under kieselguhr. B Kieselguhr alone, or on 


starch. C Kieselguhr under starch. 


each material was constant. The line A shows kiesel- 
guhr data Jf, Jd, which are independent of compres- 
sion by kieselguhr alone. The line B shows that the 
kieselguhr is subject to an apparent compressive effect 


The 


when below a starch cake, as in data Je, lg. 


precise reason for the difference between data B and 
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data A is not clear. Kieselguhr beneath a thick bed of 
kieselguhr (line 4) does not show the compressed form 
given by a similar mass of starch (line B). This may 
be due to a great difference in the lateral support 
against centrifugal stress afforded in the two packed 
materials, Starch data on line C show no change in 
packing density whatever form the laminated or 
homogeneous cake may be given. 

The linearity of each plot shows that changes in 


packing density and hence in permeability distribution 





/ 


le, 
Ly 











. é 8g 
Mass W—— 


Fig. 4. Relation between thickness of cake in small cell 


and dry weight of material in it: A Chalk alone. B Kiesel- 
C Kieselguhr under chalk. D Starch alone. 


guhr alone. 
in cakes cannot be found by thickness measurements, 
even though the mean thicknesses measured here 
were probably of a higher accuracy than the means 
obtained earlier by dyeing and sectioning cakes. 
Similar data were obtained for starch, chalk and 
kieselguhr in the small cell. In these cases the final 
sample was carefully extruded from the cell by a low 
air pressure applied to the perforated cell base. The 
sample was then sectioned for measurement of the 
layers therein. The results are shown in Fig. 4. Each 
material alone gave a straight line indicating constant 
packing density. The kieselguhr again had the pecular- 
ity of apparently closer packing when beneath a less 
permeable material. In this case the inner material 
was chalk. The results in Table 2 include data show- 
ing that the packing in the hydroextractor and the 
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small cell was the same for a given material to the 
accuracy of these measurements. Though this agree- 
ment of packing density appears to confirm the agree- 
ment found between hydroextraction permeabilities 
found in the two forms of apparatus, the confirmation 
is of minor value as the permeability tests are bound 
to be more sensitive to porosity changes than linear 
measurements of cake dimensions. The kieselguhr 
compression caused by loading with an inner layer 
of less permeable material was the same for either 
chalk or starch inner layers, although these materials 
differed considerably in permeability and hence stress 
transmitted to the kieselguhr. The apparent compres- 
sion of the kieselguhr beneath starch was not due to 
the infiltration of the starch giving an apparently 
thinner layer of kieselguhr. The layers in such cakes 
were quite easily separated, showing no macro- 
scopic mixing at the interfaces. 

The results of the permeability tests are exemplified 
in Tables 3,4 for cakes of starch and kieselguhr in 
the hydroextractor. The data are consistent in show- 
ing large increases in permeability with increasing mass 
of material in the basket, when tested at constant 
speed of revolution. This is not the expected result 
when considering that increasing the charge by adding 
another layer of cake must increase the centrifugal 
compressive stress on all layers in the initial cake. If 
the effect is only a phenomenon due to compressive 
stress in the cake, it means that the increased packing 
density and decreased permeability produced in the 
initial cake is more than offset by the higher porosity 
and permeability of the new inner layer at smaller 
radius, thus leading to an increase in mean permeabili- 
ty of the whole cake. This hypothesis requires in the 
cake a sensitivity to compressive stress, particularly 
at low stresses, which is not demonstrated by previous 
comparisons of the compressive effect in the hydro- 
extractor with loads applied directly in a compression 
cell [2], or by the insensitivity to speed of revolution 
shown in Table 3. The increase of speed from 1050 
to 1400 r.p.m., increasing the centrifugal stress by 
(1400/1050)?, produced very small decreases in the 
cake permeabilities. The above effects of radial posi- 
tion on compressive stress must influence the perme- 
abilities but it seems unlikely that the major cause of 
the large effect of charge mass can be in this field. 

A more likely cause of the phenomenon is the 
presence of a significant resistance at the cloth-cake 
interface, which is a common conjecture in filtration 


studies. Such a resistance force could arise from an 
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approximately constant effective permeability around 
the interface and thus the mean cake permeability 
would fall with increasing mass of cake, the interface 
resistance becoming less important. The dirty cloth 
on which cakes have been formed is itself of flow 
resistance only of the order 1% of the cake resist- 
ance [1], [2]. Previous tests of the effect of the 
cloth on drainage from the cake to the holes in 
the hydroextractor basket wall were consistent in 
showing that the flow through the cloth was un- 
important, either in hydroextraction or in filtration 
under directly applied load [2], though it should be 
noted that the investigation at that time was concerned 
with the source of large changes in permeability and 
the hydroextraction data were probably in error due 
to the visual detection technique used for the critical 
wetting rate [1]. The present data was considered 
far more reliable for an investigation of the cloth 
phenomenon by studying the functions of permeability 
and cake mass. 

Consider the hydroextraction permeabilities found 
by measuring the drainage time ¢ for liquid within 
a cake to move from r, tor, . The eq. (2) suggests 


the use of a plot of ¢ against log, ry r,, when using the 


same ry, 7, and r, , in order to show the effect of 


cake thickness. The plot should be linear, since C, 
would be constant, if the permeability A, were also 
constant in the equation 


(4) 
where 
ri, | 
AS 
In many tests the lines on such plots are shallow curves 
within the maximum experimental range imposed by 
the speed and geometry of the machine used. Fig. 5 
includes examples of such plots for cakes of maize 
starch, and cakes of kieselguhr. If we presume (’, to 
be defined correctly, and its value to be constant, the 
curvature of the lines in Fig. 5 is attributed to changes 
in K,. This assumption leads to the values in Table 3, 4, 
with A, increasing as W increases. If the changes in 
permeability were continuous with increasing cake 
mass, the resistance to flow residing completely in the 
cake itself, the curves in Fig. 5 should pass through 
the origin of coordinates, as shown by the dotted line 
for the region wherein no experimental data have been 
obtained. The possibility of using thin cakes in the 
hydroextractor is limited by the tendency to crack 
when supported on a flexible cloth. It appears likely 
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that the curved full lines on Fig.5 would extrapolate 
to the broken lines rather than to the dotted line. The 
former are equivalent to systems containing resistances 
to flow indicated by the intercepts on the time axis 
when the cake thickness has been reduced to zero and 
rg = 7-, the dotted lines would represent systems 
without a resistance when 7, = r, but with a rapidly 
increasing mean permeability as cake thickness in- 
creases, ¢.g. the permeability of the starch cakes at 
n’ = 1050 would have to change four-fold, mainly at 


-—--—--- -—--- — ‘ 





; 
r 10s | 


rn 
— 


c 





Fig. 5. Plot of drainage time and cake dimensions 


representing eq. (4). 


very low cake thicknesses, to provide the gradient 
the dotted 
likelihood that this rapid change will occur than that 


reduction on curves. There seems less 
the broken lines represent a condition nearer the truth. 
The gradient probably suffers an approximately con- 
stant rate of change, the curve decreasing to an inter- 
cept on the time axis. If the resistance to flow ap- 
propriate to this intercept be postulated, the question 
the cake 


thickness, e.g. due to varying compressive stress, and. 


arises whether the resistance varies with 
if so, whether its variation can be detected or analysed. 
The analysis of such a case is of perhaps small moment 
in practice but may be of value in recognising the 
cause of the phenomenon. The smooth fall in grad- 
ient as A, increases with increasing cake thickness is 
reasonable on the basis of previous knowledge of the 
dependence of permeability on compressive stress [2 |. 

Without exception the data for the present ex- 
perimental materials gave straight line plots for time ¢ 
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against either cake mass W or (7? — 7°), ¢.g. Fig. 6, 7 
Table 1, 3. 4. This 


relationship should not obtain for cakes of constant 


include the data derived from 


permeability, for which the appropriate equation is 


rh, | aXor, 


aXor?—W 


“ng 


. log, . log, (Sa) 


4° nt? kK r r 


2 
L, 
L 


where W, =2Xor2. the mass of cake material to 
fill the whole cylindrical volume of the basket of 


297 


sec / 








Fig. 6. Relation between drainage time and volume of 


hydroextractor cake [eq. (6)). 


radius ry. The broken curves on Fig.7 are hypothetical, 
being obtained from eq. (5) for selected constant perme. 
abilities. The curves A, B, C pass through the origin 
of coordinates and in contrast to the experimental 
data show an appreciable curvature. 

Also the plot of A,(W 
a straight line (Fig. 8) where “— W,” 
at ft 
the forms from Fig. 7. 


W,) against log r, 7, is 
is the intercept 


(} in Fig. 7. This is derived from eq. (4) and 


t= C,(W-W,) (6) 


where C, is a constant for a series of cakes of one 


material at one speed of revolution. Thus 


, > by 
t= C,(W+W) = 5 . log. he 
’ 
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i ke ae oe a oe a 


sec 








"400 

Moss W—e 
Fig. 7. Relation between drainage time and mass of hydro- 
extractor cake [eq. (6)]; A, B and C are hypothetical 


lines for constant permeability, obtained from eq. (5). 








Fig. 8. Plot representing eq. (7). 


whence: 


K.(W+W) = © log." (7) 


, . 
The proper relationship for constant A, is as follows 
for any value W, 


K.(W+W,) = [aX o(r2 — 732) +). slog, . (8) 


t r, 


The experimental fact of eq. (6), incorporating an 
appropriate value W,, reduces this form to eq. (7). 
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Table 5. Permeability of innermost layers in hydroextractor 


Starch: 





Layer No.* 
5 
DS ee ge ky 300 4000 500 600 TOO S00 


2-39 4-07 4:3 4:3 4-05 4-95 


Mean K for layers 2-9 is 4-6 > 


Laver A x 107 (g sec?) 


107? g/sec? 





Kieselguhr: 





Layer No.* 


3 / 


Dr Ps as Seki ab wd 200 300 400 500 
Layer K © 10 (g sec?) 11-8 51-7 42-07 44-6 


7 9 /sec? 


Mean K for layers 2-5 is 50 « 10 





* Layers built in their numerical order. 


The plot in Fig. 8 has apparently removed the time f, 
though leaving the permeability in terms of all the 
other test conditions. 

The permeability of the cakes has been analysed 
into permeabilities for the layers added in the building 
of a cake. With the usual assumptions as to flow 
conditions [1], [3] the centrifugal head developed by 
the layer of liquid from r, to rg must be balanced by 
the sum of the frictional energy requirements for the 
flow rate q appropriate to each layer of the cake 
whence this head may be expressed in terms of either 


mean or layer permeabilities as follows 


(22 n)* (72 — 19) 





Then the inner radius of the nth layer, 7,,, is the inner 


nn? 


radius r, for the cake. Whence the relationship be- 


tween the mean permeability A,, and the layer values, 


’ rn 1 
log. loge 
rs r; 
bees i) 
K, 5 


7 
loge 
- 


The permeability of a layer added to a cake was found 
from the change in mean permeability of the cake 
caused by the layer, assuming that the compressive 
effect of the thin layer was negligible. Thus eq. (10) 
was used for two layers, the initial cake providing the 
outer layer. 
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** W is the total mass of material in the cake (g). 


The data in Table 5 were obtained for series of 
tests on starch and kieselguhr in the hydroextractor. 
The initial cake layer was always of considerably lower 
permeability than the added layers. This suggests 
the presence of a high resistance in the outer layer 
of the cake. 

The data in Table 6 however for cakes of chalk and 
starch in the small cell (Fig. 2) are expressed as the 
mean hydroextraction permeabilities for a given mass 
of material in a complete cake thickness and do not 
show the low value in the first layer as given by starch 
in the hydroextractor. The data are similar for three 
methods of formation of the cake, whether the cakes 
were formed individually or in layers in the spinning 
cell, or whether they were formed by filtration in the 
cell and then spun. The three types of cake only show 


Table 6. Hydroextraction permeabilities in small cell 
10’ g/sec?, the mean value for each cake 


Chalk: 


Expressed as K 





Wig) 


Cakes formed in centrifuge 
Cakes formed by successive 
addition 


Cakes preformed by suction . 





Starch: 





W (gq) 


Cakes preformed by suction . 
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minor differences in the variation of A, with W, in- 
cluding data for A, decreasing with increasing W. 
The data in Table 6 were obtained from cakes main- 
tained under test for a considerably longer time than 
the cakes in the hydroextractor. The tests represented 
in Tabled 


a whole series of W. 


5 were of about an hour in duration for 
The tests in the small cell were 
much slower due to difficulty in cake formation and 
the production of g, # curves for each cake {3}, which 
thus entailed a total cake life of 4-6hrs. In view of 
observations in this laboratory on the variation of 
permeability with time, the data in Table 6 are consid. 
ered less reliable than those in Table 4, 5 for the present 
detailed study of permeability distribution 

The kieselguhr cakes formed by spinning in the 
small cell showed the low permeability in the first 
laver formed. The data in Tahle 7 include the mean 


permeabilities for each cake and the layer values 


Table 7. Kieselquhr in the amall cell 





Wig 2 / 6 


H ydroextraction permeability 


Shed 
val 


Mean A we 3-2 
2 13-0 


Laver*® A» 3 


Filtration perm ability 


Mean A hue 274 45 au 
Layer*® A» We 274 12-8 13-4 


12-7 





* The layer permeability is for the innermost layer, ¢.¢ 
that given for hydroextraction under mass 6¢ in the table 
(14-7 the 
changed the mean permeability 5-27 1o* for the 
lo* for the 6¢ cake 


laver, which 
ig cake 


wo*) is permeability of the 2¢ 


to the mean value 67 


deduced from the two-laver form of the equation based 
on the summation of the frictional energy requirements 
for the multi-layer cake: 


thy 
kK. (la) 


The same cakes were also tested under filtration 
conditions using the arrangement in Fig. 2 and the 
technique appropriate to eq. (3). The values obtained 
for Ay» are very near those obtained by the hydro. 
extraction equations and also show the less permeable 
first laver 

From earlier tests [2] it seemed unlikely that the 
rapid decrease of permeability with increasing radius 
a8 positions near the basket wall was due to compres. 
sive effects. As possible confirmatory evidence of this 
the porosity distribution has been measured in cakes 
of chalk and of starch. Samples were cut from speci- 


fied positions in cakes which were submerged through- 
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Radial distribution of permeability in cakes 
out their life in the hydroextractor, i.¢. the cakes were 
not spun dry before removal from the basket. The 
water content in each specimen was determined gravi- 
For 


starch the mean porosity at the front of the cake, /.¢. 


metrically and was used to derive porosities. 


at r,, was 1-035 times the mean porosity at the back 
of the cake, i+. 
a Lea and Nurse [7| apparatus for the measurement 


near rg. The dry starch packed in 





| a a a 





Fig. &. Permeability from air flow tests. 

of specific surface showed a linear relationship be- 
tween permeability to air and the porosity function 
tx (1 


milar linear function 


e)*}, as shown in Fig. 9. Assuming that a si- 
holds for the flow of 
through the packed material, the above porosity 


water 


distribution would cause a ratio of permeabilities of 
1-22: from inside to outside of the cake. This change 
of A, is insufficient to explain the high curvature 
at low cake thicknesses assumed with the dotted lines 
on Fig. 5, the latter needing the permeability distribu- 
tion of curve B on Fig. 10. 

The even gradual change of A, seems more likely 
and thus the broken lines on Fig. 5 have been accepted 
as nearer the truth. The full distribution of permeability 
with radial position is thus postulated to be that in 
curve A on Fig. 10 rather than that of curve B. The 
eurve A is based on an interface resistance defined 
as R in the centrifugal force balance 


g 


l 
3 (22n)* . (r2 — 5) ee ¥K; .q. log, y -R. (th) 


In a system with r, maintained constant, the resistanc- 
es are both functions of flow rate, g. If it be supposed 
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that the cake of thickness (7, — r,) could be maintained 


under the same centrifugal head without the resist- 
ance R, then the higher flow rate q' produced would 
be related to R by 


‘eo 


q) . log, (12) 


If the times for the flow of a volume V be 


Fig. 10. Mean permeability changes with radial position 


1 Assuming the existence of cloth cake interface resistance 


B Assuming no interface resistance 


considered asf and t! then the above force 


/ t' 
ratio is on where ty 
t ' 


lo 


(¢ -f'). This time 


difference ft, can be shown to be equal to the intercept 


on the ¢ axis in the plot of ¢ against log, r, 7, (Fig. 5) 

From eq. (1), (2) the flow rate gq when the liquid 
radius is maintained at r, is specified in terms of the 
time ¢ for the liquid laver to move from inner radius 


ry _ mm the variable ry drainage test described 


above 
l 


log, 
f 


ri) log. | ; (13) 


0 ry 
Consider a mean flow rate q,, such that the volume of 


liquid aX (rj — rj) is discharged in time ¢ 
1X(r* f ) 

aE / 
Im ' (14) 


The rate ¢ is equal to q,, when r;, is given by 


r ry 
ri )/loge|-5 
ie 
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Radial distribution of permeability in cakes 


Thus (73 — 75) is the logarithmic mean of (74 — rj 


The 


about 1-2 and thus as a good approximation 


and (r: tt ). ratio (r; rt ) (rs ry ) Was 


” » » 
(5 r’, 7 


0 L,? eo Ty) 


(72 


L 
with 7, normally near 7, 
ry aa . (17) 
The agreement between the approximate value from 
eq. (17) and the true value from eq. (15) is well within 


S31 


F 7-64 em in this work ( Table 3). both the values 


the experimental accuracy, ¢.g. with r, and 
of r, were 7-98 cm. The above relationships show that 
the time ¢ for the experimental movement of the liquid 


surface from r, tor, may be regarded as the outflow 


time for the volume between these radii for a system 


at 7, maintained at the mean of r, and r, . The ex- 
perimental drainage time ¢ may thus be considered 
as the time to be used in the relationships derived 
above from eq. (13). if the quantities Vand q therein 
be considered as the values aX "5 i) and 
aX (r7 r;_ ) t respectively. 

Thus the intercept f, on the time axis when f~=! 
in Fig. 5 is assumed to give the interface resistance as 
the fraction ¢,¢ of the total head dissipated in the 
hvdroextractor. The data for the dotted line for starch 
Fig. 5 Bon Fig. 1, 


have been reassessed by calculating the permeabilities 


at n LO50) on giving curve 


from eq.(2) appropriate to time (f — ¢,) rather than ¢. 
The curve A is obtained which shows a gradual change 
in permeability in keeping with the measured porosity 
The range of A, on 1-24 1 as 


against 42 1 for curve Band 1-22 1 from the porosity 


changes. curve A is 
distribution 

CONCLUSIONS 
No variation in the porosity of the lavers of hydroex- 
tractor cakes or the cakes in the small cell could be 
detected by measuring the thicknesses or internal and 
external radii of such cakes. The packing density was 
found to be constant throughout the cake with fine 
materials and was the same for cakes formed in the 
hvdroextractor basket or in the small cell. 

The kieselguhr layers were apparently compressed 
by either starch or chalk layers on them, while outer 
lavers of kieselguhr were not compressed by further 
added lavers of the same material. This may be due 
to differences in lateral the 


support exerted by 


materials or to the lower permeability of the starch 
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and chalk leading to higher compressive stress due yravitational acceleration (em sec?) 
to frictional drag. hydrostatic head in filtration tests (em) 
Drainage tests showed gradual decrease in perme- permeability (g sec*), A,, A,, A, fordifferent layers; 
ability with position of increasing radius in the hydro- K.. for hydroextraction; A» for filtration; K,,, for 
; on mean value 
extractor cakes. The permeability of the innermost 
; : . length of filter bed in filtration tests (em) inner 
laver was alwavs higher than that of the rest of the cake 
. . speed oft revolution (rev. sec), Ww (rev. min) 


and was nearly constant regardless of its radial position. 
' drainage rate through cake (ce sec) 


Near the outer surface, the decrease in permeability radius from basket axis (cm); r.. for basket 


was large: the value for the innermost layer of a thick surface; re for cake outer surface; r, for cake inner 


cake being four times that of the outermost laver. surface; r, for liquid inner surface: r,, 
Analy sis of the experimental results was based on for interfaces between cake lavers 
an assumed interface resistance dependent on the time (seconds) 


flow rate at the cloth-cake contact. This assessment volume (ce) 
gave a permeability variation within the cake of the mass of dry solid in cake (g) 


order of 1-24:1 from the inner to the outer laver which height of centrifuge basket (em) 

agreed with permeability distribution estimated from porosity 

. , one viscosity of liquuid (potse ) 

yravimetric determinations of the porosity. The data - ; * 

. : ‘ aa : vacking density of solids in the tested cake: (g soli 

for kieseleuhr in the small cell (Table 7) show the | 
ce cake volume) 

same permeability distribution for both filtration and 

hvdroextraction. The phenomenon is not therefore REFERENCES 

1! Borak, N. and Srorrow, J. A.; J. Soe. chem. Ind. 1950 

69 S 2) Guerra, O. Po: M. Se. Tech. thesis. Manchester 1948 


3) Harnent MoM. and Srorrow,. J. A.: Hydroextraction 


a characteristic of either hydroextraction mechanism 
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Summary The liquid vapour equilibria of two test mixtures, fs. n-decan trans-dlecalin and n-hexa 
clecane n-heptyibenzoate, were examined at low pressures; the results were checked for thermodyvnami 
consistency. It was found that the n-decane —trans-decalin system obevs Raouvts law, relative volatility 
values varving from 1-185 at 1-33 kN m l}mm Hg) to 1-36 at 101-5 KN m* (760mm Hg). The a-hexa 
decane—n-heptylbenzoate system contains an azeotrope at about 90% n-hexadecane (pressure range 
13-2400 KN m*, 1-15 mm Hg) which limits its applicability as a test mixture to evaluation of columns 
with not more than 16-18 theoretical plates 

Zusammentassung Die Flussizkeit Damptglerchgewichte der bindaren Gemische n-dekan  trans-dekalin 
und »#-hexadekan——n-heptylbenzoat, wurden bei niedrigen Drucken bestimmt; es wurde nae hygepruft 
ob die Ergebnisse thermodynamisch ino Ubereimstimmung waren. 

Es zeigte sich, da das n-dekan —trans-dekalin-System dem Raouutschen Gesetz « ntspricht und dab 
che relativen Flichtigkeitswerte von 1-185 bei 1-33 KN m? (10mm Hg) bis zu 1-36 bei LOLS KN m?* (760 mom 
Hg) variieren. Das n-hexadekan——x-heptylbenzoat-System enthalt ein Azeotrop bei ungefahr 90% + 
hexadekan (Druckgeliet 13-240 kN m 1-l5imm Hg), das die Verwendbarkeit dieses CGemisches be 


schrankt auf die Beurteilung von Kolonnen mit nicht mehr als 16-18 theoretischen Boden, 


The increasing importance of rectification at low columns interesting. Wui-iiams [19] summarized the 
pressures makes the availability of test mixtures requirements which should be met by such test 


for calibration and testing of low pressure distillation mixtures as follows 
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a) The boiling points of the components should 
be close together and should be in the range 100-200° C 
at the operating pressure; 

b) the mixtures should not deviate too much 
from ideality; 

c) any mixture of the components should be 
accurately analyzable, preferably by a simple and 
rapid method which requires only a small sample; 

d) the components should remain liquid at tem- 
peratures well below the operating temperature, 
preferably even at room temperature ; 

e) the 
chemically stable at the operating temperature; 


components should be thermally and 

f) the components should be readily available 
and easy to purify. 

However, apart from the above requirements 
relating to the nature of the components to be 
selected, a further demand as to the vapour-liquid 
equilibrium of the test mixture should be added, 
eg. the data should be accessible to control on 
thermodynamic consistency. In general, the technique 
of determination of liquid vapour equilibria is rather 
difficult, and this may be said particularly of the 
measurements of low pressure equilibria. Thermo- 
dynamic control is here a powerful expedient in 
estimating possible errors made during the deter- 
mination. 

Unfortunately, most test mixtures proposed [7], 
[11], [19] are not accessible to control on thermo- 
dynamic consistency, mostly owing to lack of ac- 
curately measured vapour pressures of the compo- 
nents. 

On the other hand, in the few cases where thermo- 
dynamic control of low pressure vapour liquid 
equilibria is possible, large systematic errors can be 
shown [3], [4], [9], [16]. 

In the present study the vapour-liquid equilibria 


The 


trans-decalin, may be used 


of two binary mixtures have been reported *. 
first mixture, n-decane 
for testing columns in the pressure range from 10 to 
760mm Hg. For pressures lower than 10mm Hg 
the n-decane—trans-decalin mixture is unsuitable in 
view of its low boiling point at these pressures. The 
second mixture, n-hexadecane—n-heptylbenzoate was 
examined in order to have a more representative 


* In the present study the dimensions of the various 
kilogram 


magnitudes are expressed in the metre second 
units system. 

In this system the unit of force is called Newton (N); 
it has the dimensions kg m/sec*. One Newton is equivalent 


to 10° dynes. 
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test mixture as regards molecular weight for vacuum 
distillation. The n-hexadecane—n-heptylbenzoate sys- 
tem may be used for calibration of columns in the 
pressure range 1-15 mm Hg. 

Both mixtures fulfil well the requirements a) to f) 
mentioned above: the components may be obtained 
by purification and/or synthesis from commercially 
available base materials; analysis may be performed 
quickly and accurately by means of refractive index 
measurement; the components are liquid at room 
temperature and stable at distillation temperatures. 
The n-hexadecane—n-heptylbenzoate mixture, how- 
ever, contains an azeotrope and deviates rather much 
from ideality; not more than 16-18 theoretical plates 
can be accurately measured with it. 

In the second paper in this series, the results 
of testing a Vigreux column and a rotating-band 
column with the binary mixtures mentioned over a 
wide range of pressures will be described. 


CHeck OF THERMODYNAMIC CONSISTENCY 

In recent literature on the checking of vapour-liquid 
equilibria for thermodynamic consistency several 
methods of correlation have been given [5], [10], [18], 
[20], frequently based on special integrations of the 
Duhem-equation. A more general method, which is 
very simple in application too, is described by Rep- 
Lich and Kuster [12] and by Covntson and He. 
RINGTON [6]. 

The method consists in calculating the ratio of 
activity coefficients from the «-values and vapour 
pressure ratio by 


(1) 


after which the log y,/yg values obtained are plotted 
against the mol fraction of volatile component in the 
liquid phase. In the case of thermodynamic consistent 
results the following condition should be met: 
1 

| log " dx=0, (2) 

. 42 

0 
so the net area between the curve obtained (‘‘ex- 
change function’’) and the abscissa should be zero. 

For accurate work, activity coefficients should be 

corrected for gas law deviations and pressure effects 
on the liquid. These corrections may be incorporated 
in the vapour pressures of the pure components 
according to Benepict et al. [2} by: 


P, “ 
Pt ~ 2.303 RT & 


(V, a B,) (P— P)) | . (3) 


log 
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The second virial coefficient B in this equation is 
calculated from critical data by BerTHeLot’s equa- 


tion: 
9 T, 
B RM |) _¢ 


Is P 


(4) 


T.\2 
- 


Vapour Ligurp EQuILIBRia 


materials 


Base 
n-Decane was obtained by dehydratation of a commer- 
cial decanol-1 and subsequent hydrogenation of the 
The base material, the olefinic 


resultant n-decene. 


intermediate product and the final product were 
distilled in a 60-plate Oldershaw column. Only the 
fractions differing less than 0-1° C in Cottrell boiling 
point were combined. 

Trans-decalin was obtained from commercial 
decalin by fractional distillation. 

In order to make the yield as high as possible, 
the cis-decalin present was first converted to ftrans- 


AICI, 


according to Sever and Yur [15]. 


decalin by isomerisation with about 20% 


catalyst at 50°C 
The reaction product was sharply fractionated in a 
60-plate Oldershaw 


column; again, only the frae- 


less than 0-1°>C in Cottrell boiling 


tions differing 
point were combined. 

n-Hexadecane was a commercial product, purified 
from a mixture of 


by repeated recrystallization 


benzene-methylethyl ketone: the recrystallized pro- 

duct was then distilled under reduced pressure 
n-Heptylbenzoate was synthesised by heating an 

with both 


excess of n-heptanol benzovichloride 


materials were previously purified by distillation. 


The 


rectification in a 


reaction product was subjected to vacuum 


Vigreux column (equivalent to 


about 15 theoretical plates). Only the fractions dif- 


fering less than ©-0001 in refractive index were 


combined. 
Physical properties of the materials obtained 


are given in Table |. The large difference in refractive 


Table 1. 


Physical properties of materials used 





Boiling point Melting point 


Pal) 
” D 
( 


1-4118 
1-41 189 
1-4696 
1-46968 
1-4692 
1-4347 
1-43448 
1-4917 


n-decane 
Litt. [1]. 

t-decalin 
Litt. [14] . 
Litt. [16] . 

n-hexadecane 
Litt. [1]. 


n-heptylbenz vate 
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index between the components of each of the two 
test mixtures offers the possibility of analysis by 
this magnitude. To this end refractive index deter- 
minations of a number of mixtures of the components 
were made; the data obtained are given in Tables 2 


and 3. 
Table 2. Relation between refractive index and composition in 


the system n-decane —t-decalin 





Mol. fr. 


n-le cane 


Mol. fr. 


nao 
) 
u n-cde cane 


"Dp 


0-284 
0-083 


1-000 1-4510 
1-4639 


1-4696 


0-607 


0-436 0-000 





3. Relation between refractive index and composition in 


the system n-hexadecane — n-he pt ylbenzoate 





Mol. fr. 
n. Cr 


” Mol. ir. m” 
"D 0 ” D 
ms 
204 


1-4714 
“SIS 1-4767 
TO 1-4826 
592 1-4017 

0-304 


1-000 


0-125 


O-O00) 





Apparatus 
Most of the vapour-liquid equilibrium determinations 
were performed with the Othmer-type still sketched 
in Fig. 1. 


or overheating of the vapours, the heater for 


To avoid systematic errors by condensa- 
tion 
compensation of heat losses was assembled around a 
jacket through which the vapours from the equi- 


librium chamber were passed off. 


In order to avoid bumping of the liquid at lower 
pressures, the liquid chamber was filled with Ni- 


chromium Fenske helices. 


This apparatus is, however, not suitable for use 
at pressures of about | mm Hg as then the pressure 
drop in the vapour line is too large. Therefore an 
apparatus as described by WILLIAMS [19] was used for 
this pressure range. It proved necessary to modify 
the apparatus slightly as it appeared from meas- 
urement of the boiling point of »-hexadecane at 
1mm Hg that there was still a considerable pressure 
gradient in the vapour line. The upper part was 
arranged like a reflux condenser, with a separate 


pressure exchange tube to the receiver part (Fig. 2). 
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Vapour pressure determinations of n-decane and 
trans-decalin were made in a Cottrell boiling point 
apparatus (liquid capacity about 20 ml) (Table 4); 


! TO MANOSTAT 
THERMOWELL 


ty 
| 


HEATING | 
MGs 


8 


INSULATION 









































‘MEATERWELL 


Fig. 1. Othmer Still. 


Table 4. Vapour pressures of n-decane and trans-decalin 





n-decane trans-decalin 


P t - P t Pp 
kN /m? C kNim kN in? C kNim?* 


179-1 
175-6 76-25 
171-9 69-6 
173-5 72-2 
177-1 78-9 
180-6 85-5 
183-7 92-2 
186-7 98.9 
188-4 102-9 


129-7 82-9 
138-3 


143-0 


60-1 1-41 
68-55 2-10 
80-9 3°56 
139-0 91-75 5-45 
132-3 100-0 745 
125-7 106-45 945 
119-0 114-4 12-60 
112-3 87-6 4-64 
105-7 188-4 102-9 

102-9 185-5 96-3 

182-4 89-5 


184-1 
186-9 
188-4 
187°] 
185-0 
182-9 
180-6 
178-2 
175-8 
174-7 


56-7 1-41 
64-65 2-10 
76-15 3-56 
86-25 5-45 
93-9 7-45 
99-85 9-45 
107-15 = 12-60 
82-4 4-64 
174-7 102-9 
179-7 = 116-4 
181-9 123-0 
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the data on n-hexadecane and n-heptylbenzoate 
were obtained by distilling the pure materials in 
the equilibrium still of Fig. 2 at various pressures 
(Table 7). 

At the lowest pressures the pressure in the appara- 
tus was kept constant by a manostat as described 
by WrtiiaMs [19]; at higher pressures a pressure 
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1] 
Fig. 2. Othmer Still, used at about 0-13 kKN/m*. 


regulator of the design given by Topp [17] was 
used, 


n-Decane—trans-decalin binary system 


The data obtained at atmospheric pressure with the 
equilibrium still of Fig. 1 are collected in Table 5. 
In checking these data for thermodynamic con- 
sistency, the following values were used in the cal- 
culation of the correction for gas law deviations ete. 
according to eq. (3) and (4). 





Molal vol, 


m®/k mol 


Crit. temp. Crit. press. 
°C MN/m?* 


m-decane ...... 344 2-09 0-234 
t-decalin ...... 2- 0-183 
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atmospheric pressure 


trans-decalin 
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pressures between pure n-decane and two mixtures 


of n-decane and trans-decalin. These measurements 





Composition, mol. 


fraction n-decane log a log * 
Liquid Vapour "s 
0-241 0-305 0-140 0-134 
0-354 0-432 0-143 0-134 
0-351 0-427 0-140 0-134 
0-338 0-405 0-130 0-134 
0-312 O-387 0-143 0-134 
O-31L5 382 0-127 0-134 
0-305 O-370 0-127 0-134 
0-300 O-370 0-137 0-134 
362 0-441 0-137 0-134 
3Ol 0-460 0-124 0-134 
0-445 0-527 0-148 0-134 
O65 0-641 O-137 0-133 
O-612 676 O-121 0-133 
0-604 0-74 127 0-133 
O-7T17 776 0-137 0-133 
OSLO O-S55 140 0-133 
O-S2D5 O-s62 124 0-133 
0-691 O70 0-127 133 
0-065 O-087 0-140 0-134 
0-218 276 137 134 
O-199 245 114 134 
O-195 O-245 0-127 0-134 
0-165 210 130 0-134 
105 O-137 134 134 


0-006 
0-009 
0-006 
0-004 
0-009 
0-007 
0-007 
0-003 
0-008 
O10 
0-000 
0-004 
o-O1l2 
0-006 
O-004 
0-007 
0-009 
0-006 
0-006 
0-003 
0-020 
0-007 
0-004 
0-000 


were carried out in two Cottrell boiling point appara- 























tuses already mentioned above. Temperature dif- 
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Fig. 3. Exchange function of n-decane—trans-decalin at 


atmospheric pressure, 


ferences were determined with the aid of a differential 
thermocouple and a sensitive galvanometer (enabling 


readings within 0-01-0-02° C). The results obtained 





By applying this correction the 
log of vapour pressure ratio of the 
the 
data in Table 4 (ranging from 0-140 


components calculated from 


to 0-143) is reduced by about 5°. 

The “log of the corrected vapour 
pressure ratio and the “log of the 
ratio of activity coefficients cal- 
culated are given in Table 5 too. 

From the plot of the atmospheric 
data in Fig.3 it may be concluded 
that the results obtained are thermo- 


dynamically consistent; the system 





behaves ideally at least within the 
limits of experimental error. 

As the vapour-liquid equilibrium 
at atmospheric pressure proved to 











are collected in Table 6 and shown in Fig. 4. 





T ——— ] - <7 - T 








be ideal, only a few check deter- 07 | | | | | | | | | | | | | | 
minations were carried out at re- o 7 2 J * 5 6 7 8 9 0 Mf 2 19 % 
luced -ssures, viz. by measuring PRESSURE kN/m2 
Gucea pressures, Utz. 0) = Fig. 4. Boiling point differences between n-decane and mixtures of n-decane 


boiling point differences at various 


and trans-decalin. 
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Table 6. Temperature intervals between boiling points of This test [13] which is highly sensitive to diverg- 
n-decane and n-decane —trans-decalin mixtures ence from the ideal, shows that the differences 


0 mol % n-decane 50-0mol %o n-decane 93-8 mol % n-decane between calculated and experimental total pressure 





: are very small and within the erro xperimenté 
kNim? #,°C  kNim® §¢,°C | kNim®~—s 6, °C A reer of experimental 
observation (about 0-02 kN/m?). 


1-49 The system n-decane—trans-decalin being ideal, 
2.8] 9. “ 2 plate numbers of columns distilling under total 
4-45 2.36 Bb 26 reflux may be calculated with Fenske’s formula [8]. 
86 2-7 7-36 . As the x-value does not vary much with concentration, 
“30 3- 9-55 an average value may be used at every specified 


20 . 2. . pressure, as follows: 





5-19 2. , 2. 2 Pressure 
a-value 





kN/m? mm Hq 

From the data with the mixture containing 
93-8 mol®% of n-decane, dtdx values at x = 1 in 10 
the ¢— x diagram may be calculated. Repticn and 9.67 »”) 
Kuster [12] have shown how to calculate « values 6-67 50 
from these data. Taking the values from Fig. 4 13-33 100 
at resp. 2-67, 6-66 and 13:33 kN m? (20, 50 and LOLS 760 


100mm Hg), the following results are obtained: 








Recently, Struck and Krxney [16] and Fenske, 
Te. 4s 6 « 2-67 13-33 Myers and QvuicGLe [7] published data on the 
mag. 2 ee es : 100 vapour-liquid equilibrium n-decane—trans-decalin over 
dijdx (Fig.4)... . 3 5 58 a wide range of pressures. They found z-values 
d log P, dt (derived 
from Table 4). . . 0-0202 0-O178 0-0162 
ee De ah ee 1-21 1-26 1-28 


differing greatly from the values presently reported, 
viz. at LOmm Hg, « = 1-11, at lOOmm Hg, « = 1-22, at 
760mm Hg, «= 1-30. The authors did not test their 


ideal « from vapour 


- . ° : data for thermodynamic consistency, however, so that 
pressures*. ... . 1-205 1-25 1-30 Bes : . , . 
large systematic errors in the data given remained 

* Corrections for divergence from ideal gas laws may be — ynobserved. Applying the criterion of Repiicn ete. 
ignored in view of low total pressure. Pe 





to some of their data, the result presented in Fig. 5 


The above calculations show that the system is obtained. It is seen that the log y,/y,-function 


n-decane—trans-decalin behaves ideally at lower press- can be represented by a straight line, running parallel 
oi i had 
/Pe 0-02. 


This can also be demonstrated by comparing the This shows that systematic errors are present in the 


ures too. to the abcissa at a distance of about log y 


measured total pressure of the mixture containing determinations made; the result, however, suggests 
50 mol % n-decane with the calculated pressure too that the n-decane—trans-decalin system behaves 
according to Raouts law. ideally. This may also be shown by calculating 

The calculations yielded the following results: total pressures according to Raovuut’s law from 





At for mixture Boiling pt. Sinein cient 
Boiling pt. 50-0 mol %o mixture pe F Calculated Difference 
Le 50-0 mol pressure calc—obs. 


(Fig. 4) n. Cy Cc tr.-dec. 
kN in mm Hg / } Cc oN kN/m? KN/m? 


Pressure observed 
n. Cro n. € 


2-67 2 69-8 71-75 2-02 2-4: 2-67 0-00 
6-67 . 91-2 93-90 . 6-68 + 0-01 
10-66 102-9 106-0 2- 10-67 +0-01 
13-33 108-6 111-9 5. 13-31 — 0-02 
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boiling points and compositions of the binary systems 
given by Srruck and Kiyney. As appears from 
calculated and experimental 


the ratio between 


pressures—also represented in Fig. 5 as a function 
of liquid composition—the system is actually almost 
ideal. The deviations found between actual and 


calculated pressures are invariably smaller than the 


PRESSURE 
SYMBOL mm hg AUTHOR 
° 760 FENSKE ET Al 
x 50 STRUCK ET Al 
° 20 ” 
P CALCULATED 10 F 
P EXPERIMENTAL 
y T 7 7 
102) ’ + , , > > - + +- — 








4 








0? 02 0) 06 05 06 07 08 09 10 


X. MOL FRACTION n. DECANE IN LIQUID 


Fig. 5. Correlation of literature data on the n-decane 
trans-decalin equilibrium ([{9], [16)) 


limits within which the pressure on the system was 
kept constant during the determinations (about 
0-5 mm Hg). 


n-He xadecane —n-he ptylbenzoate binary system 


The vapour-liquid equilibrium of the n-hexadecane 
n-heptylbenzoate system was measured at three 
pressures: 0-13, 0-49 and 24)kNm* (1, 3-7 and 
15mm Hg, respectively). For determinations at 
0-49 and 2-0 kKN’m* the Othmer apparatus of Fig. | 
was used; the measurements at 0-13 kN m* were 
carried out in the equilibrium still of Fig.2. The 
data on the vapour pressures of the components 
and the vapour-liquid relations obtained are collected 
in Tables 7-10. 


Chemical 
Engineering Science 


Table 7. Vapour pressures of n-hexadecane and n-heptylbenzoate 





n-he radecane n-hept ylbenzoate 


t P t P t P t P 
Cc Nim? Cc Nm Cc N/im?@ Cc Nim? 


111-1 I8o))—-130-2 550 61108) «6113 107-2 96 
109-8 8163) 125-9 430 «101-2 59 «(106-3 92 
7-8 144) 134-9 730 #1156 162 1048 82 
105-7 128s 615 112-1 133. 103-0 74 
116-0 256 137-8 860 131-7 395 123-9 266 
13-1 208 140-1 905 1281 318 139-5 640 
110-7) «= 179s«d136-9 800 «121-90 236 «148-0 933 
1Ww7-0 139 i520) «1640)—CO- 180 156-8 1440 
103-5 115 158-7 280 127-8 320 170-5 2390 
100-6 92 1904 7250 1223 295i 176-5 2970 
OS-S so Is 2oso 142-1 730 182-0 3670 
07-4 73 1490S 1450) 135-9 545 
31-5 0« «65854200 BB 400 
122-3 368 170-9 §=3620 = «108-0 68 
120-1 = 328 101-0 62 
133-9 695 125-8 420 111-4 27 








Table 8. Vapour-liquid equilibrium n-hexadecane —n-heptyl- 
benzoate. Experimental data at 133 N/m* 
t ‘om position 
mol. fraction Tempe- P y 
n-hezadecant rature log « log Pp. log ie 
2 42 
Liquid Vapour Cc? 
(822 0-828 105-1 O-O1S 0-181 0-163 
760 782 105-1 0-054 0-181 0-127 
719 0-744 105-1 0-055 181 0-126 
712 0-742 105-1 0-065 181 0-116 
0-640 0-686 105-3 0-090 0-180 0-090 
0-500 0-603 105-9 0-182 0-180 + 0-002 
477 0-587 106-0 0-192 0-179 0-013 
417 0-533 106-4 0-203 0-179 0-024 
0-230 0-403 108-2 0-354 0-177 O-177 
212 0-391 108-5 O-377 O-177 0-200 
0-186 0-351 108-8 374 0-176 0-198 
170 0-320 100-1 0-361 0-176 O-1S85 
0-163 0-302 109-2 0-346 0-176 0-170 
0-135 0-268 109-6 0-370 0-175 0-195 
0-082 0-209 110-6 0-471 0-174 0-297 
0-096 0-188 110-3 0-340 0-174 0-166 
0-086 0-183 110-5 0-376 0-174 0-202 
0-073 0-153 110-8 0-360 0-174 0-186 
0-057 0-127 11l-l 0-382 0-174 0-208 
0-794 0-808 105-1 0-038 O-181 0-143 
0-812 0-813 105-1 0-003 0-181 0-178 
0-842 0-850 105-1 0-026 0-181 0-155 





* Smoothed from experimental data. 
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Table 8. (Continued) Table 9. (Continued) 





Composition Composition 

mol. fraction Tempe- P mol. fraction Tempe- P, 

n-he radecane rature log log P n-hexadecane rature log « log P 
2 Ye 2 


Liquid Vapour °os Liquid Vapour 


0-888 0-886 105-2 0-009 0-180 . 0-140 0-253 32: 0-318 0-157 0-161 
0-925 0-925 105-4 0-000 0-180 00S «O10 sale 0-298 0-156 0-142 
0-847 | 0-847 105-1 0-000 | 0-181 ; 0-038 0-083 34: 0-360 0-156 0-204 
0-847 = 0-847 105°1 0-000 = O-181 — | = ed oc | Om O-183 
0-784 | 0-790 105-1 0-014 0-181 oc | OCs | 0123 | o150 0-036 
0-767 0-779 = 105-1 0-029 0-181 0-152 0-631 0-694 28- 0-123 0-159 0-036 
0-909 = 0-896 105-3 0-065 0-180 0-245 oon | on _ orm | om O08 
0-886 0-900 105-2 0-064 0-181 0-171 owe | Or pay oces | O18 one 
v-900 | 0-900 105- 0-750 0-778 27-5 0-067 0-160 0-093 
0-903 0-900 105: 0-798 0-807 27: 0-025 0-160 0-135 


0-000 0-181 0-181 
0-015 0-181 0-196 
0312 0-491 107: 0-342 860-846 27: 0-012 0-160 0-148 
0-271 0-439 107- 
0-586 0-662 105- 0-140 0-180 0-040 
0-133 0-293 109-6 0-431 0-175 0-256 0-973 0-973 7-6 0-000 = 0-160 0-160 
0-073 = O-189 110-8 0-472 0-174 (0-298 0-968 0-972 27-6 0-059 =—-0-160 0-101 
pi po 110 0-431 0-174 9.en7 0-936 0-936 97. 0-000 0-160 0-160 
0-068 0-152 110-9 0-390 0-174 0 0-134 = 0-228 32-6 0-281 0-157 0-124 


0-057 0-133 111-1 0-404 0-174 0 


0-328 0-178 0-150 

7" 7 . \o -16 .).59 
0-323 0-177 0-146 0-873 0-871 0-008 0-160 0-152 
0-896 0-896 27- 0-000 0-160 0-160 


. 
. 
0327 | 0492 | 1071 «= 0-300 o-178. | «So ogg «82 828127 oe) Sa | eee 
2 
7 
5 


57 
16 


30 





“2 
2 
2 
_ Table 10. Vapour-liquid equilibrium n-hexadecane — n-heptyl- 
benzoate. Experimental data at 2-0 kN/m* 








Table 9. Vapour-liquid equilibrium n-hexadecane —n-heptyl- Composition 


benzoate. Experimental data at 490 N/m?* mol. fraction Temp. oe 
n-hexradecane P. . 


42 





Composition 
mol. fraction Temp P, Liquid Vapour 
> 


n-hexadecane log a log I - 

3 0-076 0-140 164-6 0-297 0-146 0-151 
0-116 0-197 163-9 0-272 0-146 0-126 
O-181 0-287 163-0 0-260 0-146 0-114 
0-196 0-301 162-8 0-248 0-146 0-102 
0-205 0-322 162-7 0-267 0-146 0-121 
O14 | O284 | 1336 0-207 | 0-157 O10 9-186 0-293 162-9 0-258 0-146 0-112 
0-137 0-214 132-6 0-233 0-157 0-076 


Liquid Vapour c® 


0-136 ()-225 132-6 0-267 0-157 0-110 
0-136 0-229 132-6 0-277 0-157 0-120 


0-282 0-383 161-5 0-199 0-146 0-053 
0-132 0-238 132-6 0-312 0157 0-155 0-385 0-480 160-4 0-170 0-146 0-024 
0-133 0-245 132-6 0-326 0-157 0-169 
0-167 0-275 132-2 0-277 0-157 0-120 
0-176 0-305 132-1 0-312 0-157 0-155 
0-244 0-373 131-2 0-267 0-158 0-109 
0-250 0-377 131-2 0-258 0-158 0-100 
0-303 0-411 130-5 0-207 0-158 0-049 


0-440 0-533 159-8 0-161 0-147 0-014 
0-470 0-556 159-5 0-149 0-147 0-002 
0-486 0-577 159-3 0-158 0-147 O-O11 
0-404 0-493 160-2 0-158 0-146 0-012 
0-565 0633 158-4 0-124 = 0-147 0-023 
0418 0-491 160-0 0-130 0-146 0-016 
0-434 0-526 129-3 0-161 0-159 0-002 0-383 0-484 160-4 0-179 0-146 0-013 
0-424 | 0-514 129-4 0-158 | 0-159 0-001 0-222 0-335: 162-4 0-246 «©0146 0-100 
0-301 | 0-416 = 130-5 0-220 | 0-158 0-062 0-095 0163 1643 0-267 0-146 0-121 
0-200 0-416 130-7 0-241 0-158 0-083 0-886 0-888 156-5 0-008 «0-148 0-140 
0-418 0-518 129-5 0-176 0-159 0-O17 0-800 0-824 156-5 0-040 0-148 0-108 
0-335 0-450 130-3 0-212 0-158 0-054 0-760 0-784 156-7 0-059 0-148 0-089 
oss | OS (1206 0238 | 0168 0073 9.619 0-680 158-0 O-116 0-147 0-031 
O341 | 0461 = 130-2 0-218 | 0-158 0-060 0-565 0-638 158-4 0-132 0-147 0-015 





* Smoothed from experimental data. * Smoothed from experimental data. 
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The ratios of the activity coefficients calculated carried out on a mixture with about 86-5 mol % 


from these experimental data are plotted in Fig. 6-8. n-hexadecane in a Vigreux column (about 10 theoret- 


It is seen that there is considerable spreading in the _ ical plates). 
I P 


results obtained. In drawing the curves through Under total reflux the top composition was 


the points the position of the azeotrope in the system determined as a function of the pressure (Table 11). 
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Fig. 6. Exchange function of n-hexadecane —n-heptyl- Fig. 7. Exchange function of n-hexadecane —n-heptyl- 
benzoate at 0-13 kN/m?. benzoate at 0-49 kN/m!?. 


4 


should be exactly known. This cannot accurately be Extrapolation and interpolation of these data gives 


determined from the equilibrium measurements them- — the following azeotropic compositions: 





selves, however; therefore distillation tests were 
Pressure Azeotropic comp. 


Table 11. Determination of az otropic composition in the kN/mé@ mol % n. Cy, 


system n-hexadecane —n-heptylbenzoate by means of distilla- 
0-13 
; - 0-49 
Pressure kN/m? Composition mol, “o n. Cy 
20 


tion in a 10-plate Vigreux column 








Top Bottom Top Bottom 


The “exchange function’’ was then laid through 
these azeotropic points so that the area condition 
is complied with. It appears that the experimental 
values, especially at low concentrations of n-hexa- 


1-45 1-60 80-1 
O-57 1-00 
0-40 0-80 


0-80 O85 
1-35 1-45 86-6 decane, exhibit relatively great deviations from the 


1-05 1-20 5 curve drawn; if, however, the nature of the meas- 
.o* “AS 5-5 . . 
Os oe < urements is taken into account, (low pressures, 
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irregular boiling), the general result is not unsatis- 
factory. 

In order to calculate plate-numbers at total 
reflux, smoothed values of « were derived from the 
lines drawn in Fig.6-8 and the vapour pressure 
the components. 
in Table 12, 


ratio of From these «-values, 


recorded curves were calculated from 


% 
L005, 
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Fig. 8. Exchange function of n-hexadecane 


benzoate at 2-0 kKN/m?*. 


Table 12. 


n-heptylbenzoate smoothed values at various pressures 


Vapour-liquid equilibrium data of n-hexadecane 





O13 kN ni* O49 kN im? 20kN/mé? 


7 4 x y y 


0-000 
0-195 
0-324 
0-423 
0-505 
0-582 
0-656 
0-730 
0-810 
0-899 
1-000 


0-000 
0-209 
0-342 
0-443 
0-523 
47 0-595 
1-31 0-662 
1-19 0-735 
1-07 0-810 
0-975 0-898 
O-89 1-000 


0-000 
0-186 
0-313 
0-413 
0-498 
0-579 
0-655 
0-730 
0-810 
0-899 
1-000 


—wem 8 tS tS 


0-99 


0-92 0-93 
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which the number of plates between top and bottom 
(differential plate number) can be directly derived 
from the 
products found (Fig. 9). 


refractive indices of top and bottom 
This figure shows that, 
owing to the existence of the azeotropic point, not 
more than 16-18 theoretical plates can be measured 
with the n-hexadecane 


20 
"D 


n-heptylbenzoate system. 


+ 


049 kN/m2 


+ + —+ 


0133 kN/m2 


+444 








144 . 
0 2 4 16 18 AN 
Fig. 9. Refractive index VS differential plates in the 
system n-hexadecane 





n-heptylbenzoate at various 


pressures. 


NOTATION 


second virial coefficient m*/k mol 
(see footnote on page 165) 

N/m? 

N/m? 

N/m? 

Nm/k mol °K 


abs. temperature K 


total pressure 

corrected vapour pressures [eq. (3) | 

vapour pressure 

gas constant 

molal volume m*/k mol 
mol fraction volatile component in liquid 
mol fraction volatile component in vapour 
relative volatility 
activity coefficient 

Indices: 

! volatile component 


2 — heavy component 
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Vacuum distillation I 
Performance of a Vigreux column and a spinning band column over a wide range of pressures 
F. J. ZUIDERWEG 
Koninklijke Shell-Laboratorium, Amsterdam 


( Received 6 Se ptember 1951) 


Summary—Performance figures referring to a 11 mm diameter Vigreux column and a 6 mm diameter 
spinning band column over a pressure range of about 0-13 KN’m*-101-5 kN/m* (about 1-760 mm Hg) 
are presented. 

In order to gain an insight into the mechanisms of material transfer involved, analyses of separate 
vapour and liquid phase resistances were made. 

In the case of the Vigreux column it was found that liquid phase resistance was about half the total 
resistance at atmospheric pressure. A reduction in pressure causes the liquid phase resistance to be in- 
creased by increasing liquid viscosity on the one hand and to be decreased by the higher vapour rates on the 
other hand. These two compensating factors cause a total efficiency loss of the Vigreux column of about 
25% when the pressure is reduced from atmospheric value to about 0-13 kN/m? (1 mm Hg). At this top 
pressure the maximum pressure drop along the Vigreux section is about 0-3 kKN/m* (3-4 mm Hg) per metre. 

In a spinning band column the main resistance to material transfer is offered by the vapour phase. 
This resistance is reduced by creating turbulence through rotation of the spinning band. The present 
study disclosed the existance of two limiting values of rotational speed: a lower limit below which rotation 
does not yet cause turbulence, and a higher limit above which back mixing in the liquid phase causes the 
separating power to be independent of—or to decrease with—speed of rotation, 

Due to the fact that the extent to which turbulence can be effected in the vapour phase proves to 
depend on rotational speed and vapour density to the same degree, the sensitiveness to pressure variations 
of the spinning band column is rather great. Reduction of the pressure from atmospheric value to about 
1-3 kKN/m? (10 mm Hg) causes an efficiency loss of about 70% ; in fact at this pressure the separating power 
of the spinning band column was not better than that of a wetted wall column of the same dimensions. 


Zusammenfassung Es werden die Leistungszahlen einer Vigreux-Kolonne, Durchmesser 11 mm, und einer 
Drehbandkolonne, Durchmesser 6mm, Uber einen Druckbereich von ungefahr 0,13 kKN/m* —101,5 kKN/m?* 
(1—760 mm Hg) gegeben. 

Um eine Einsicht in den Mechanismus des Stoffaustausches zu bekommen, wurden die Widerstande 
der Dampf- und der Flissigkeitsphase separat analysiert. 

In der Vigreux-Kolonne ergab sich der Widerstand der Flissigkeitsphase zu ungefihr der Halfte 
des totalen Widerstandes bei atmospharischem Druck. Druckverminderung veranlaBt einerseits eine 
Zunahme des Widerstandes der Flissigkeitsphase infolge erhéhter Flissigkeitsviskositat und andererseits 
eine Abnahme infolge héherer Dampfgeschwindigkeiten. Diese zwei einander kompensierenden Faktoren 
verursachen einen totalen Leistungsverlust der Kolonne von ungefahr 25%, wenn der Druck vom atmo- 
spharischen Wert bis zu ungefahr 0,13 kN m* (1 mm Hg) erniedrigt wird. Bei diesem Druck am Kopf der 
Kolonne ist der maximale Druckverlust entlang der Kolonne ungefahr 0,3 kKN/m? (3——-4 mm Hg) pro Meter. 

In einer Drehbandkolonne wird der Hauptwiderstand gegen den Stoffaustausch durch die Dampf- 
phase geleistet. Dieser Widerstand wird mittels einer durch Drehung des Drehbandes verursachten Wirbel.- 
bewegung herabgesetzt. Die Untersuchung ergab die Existenz von zwei beschrankenden Umdrehungs- 
geschwindigkeitswerten : einen niedrigen Grenzwert unter welchem die Drehung noch keine Wirbelbewegung 
verursacht und einen héheren Grenzwert uber welchem Rickmischung in die Flissigkeitsphase die Ursache 
ist, daB die Trennungsfahigkeit von der Umdrehungsgeschwindigkeit unabhangig ist oder damit abnimmt. 

Da das MaB, in welchem eine Wirbelbewegung in der Dampfphase erreicht werden kann, im gleichen 
MaBe von der Umdrehungsgeschwindigkeit und Dampfdichte abhangig ist, ist die Druckanderungs- 
empfindlichkeit der Drehbandkolonne ziemlich gro8. Herabsetzung des Druckes von atmospharischem 
Wert bis zu ungefahr 1,3 kN/m* (10 mm Hg) verursacht einen Nutzeffektverlust von ungefahr 70% ; 
tatsichlich ist die Trennungsfahigkeit der Drehbandkolonne bei diesem Druck nicht besser als die eines 
leeren zylindrischen Rohrs von gleichen Abmessungen. 
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INTRODUCTION 

Though fractionation at subatmospheric pressures 
may be employed as a means of obtaining a sharper 
separation (e.g. by increasing the relative volatility of 
the components), vacuum distillation is generally used 
in order to lower the boiling points of the components 
when heat-sensitive materials are to be distilled. The 
maximum effect to be obtained from the decrease in 
pressure applied depends, however, upon the pressure 
drop occurring between top and bottom of the frac- 
tionating column used. This largely governs the 
selection of types of fractionating column in vacuum 
distillation. 

When the pressure in the reboiler should be below 
13-3 kN m?* (100mm Hg), bubble plate columns cannot 
well be used unless the number of plates is very low; 
some type of packed column must then be resorted to. 
When the pressure in the reboiler must be still lower 
(e.g. below 1-33 kN m? [10 mm Hg}), the pressure drop 
of this type of column is too large too [28], so that in 
this pressure region use should be made of open tube 
columns. 

If the thermal hazard of the material distilled is 
still too high when using these columns, molecular or 
semi-molecular distillation may be applied; the pro- 
cedure to obtain a sharp separation with this type of 
apparatus is, however, not simple. 

Thus far, published information on the perform. 
ance of low-pressure fractionating columns has been 
scarce. Moreover, practically no efforts have been 
made to give an interpretation of the pressure effects 
found. Perhaps the best study in this respect was 
made by ScHAFFNER, BowMaAN and Cou. [26] on 
their empty tube column with alternating hot and 
cold sections. 

Discussion of published results is further hampered 
by the fact that sometimes separating power data have 
been calculated with thermo-dynamically inconsistent 
equilibrium data [2], [28]. 

In general it may be concluded, however, that 


packed columns lose distinctly in separating efficiency 


when the pressure is reduced. FELDMAN [16] et al. 
showed a Podbielniak Heligrid column to have a 50% 
lower separating power at 6-7 kN m? (50 mm Hg) than 


at atmospheric pressure. 


* In the present study the dimensions of the various 


magnitudes are expressed in the metre—kilogram—second 
units system. 

In this system the unii of force is called Newton (N); 
it has the dimensions kgm sec*. One Newton is equivalent 


to 10° dynes. 


On the other hand, open tube columns seem to 
retain their separating power fairly well. Measure- 
ments by WILLIAMS [34] at atmospheric pressure and 
at 0-13 kN m?* (1 mm Hg) show only a small efficiency 
loss. Bretwiscu [4] determined the separating effect 
of a number of wetted wall columns of different dia- 
meters under laminar flow in the vapour phase. At 
pressures of 13kN m*—4-0 kN m? the results cor- 
related well with the theoretical separating power, as 
is the case at atmospheric pressure [23], [32]. 

Sometimes columns with rotating elements are re- 
commended for use in vacuum fractionation. Bircn, 
NaTHAN and Gripp [3] describe a rotating band co- 
lumn of large diameter and give test data at atmo- 
spheric pressure. Low-pressure separating power data 
were not determined; distillation curves obtained at 
low pressures, however, suggest a reasonable separat- 
ing power. 

In this study the performance of two types of open 
tube column frequently used in vacuum distillation 
was investigated. Performance figures for a Vigreux 
column and a rotating band column over a wide range 
of pressure are given. They were determined with the 
low-pressure test mixtures n.decane-transdecalin and 
n.hexadecane-n.heptylbenzoate described in the com- 
panion paper [36]. 

An interpretation of the pressure effects found has 
been given. In order to facilitate the discussion of 
the experimental results, a comprehensive survey of 
the principles that govern material transfer in the 
columns tested will be given in the following para- 
graph. 

THEORETICAL 

In the evaluation of laboratory distillation apparatus 
the separating power is generally expressed in the 
theoretical plate number. The influence of operating 
variables on the overall performance of the apparatus 
is thus easily recognized and, moreover, the plate 
number found may be used in simple approximate 
calculations to investigate the utility of the apparatus 
in actual distillations. In the present work the same 
criterion will therefore be applied too. 

In the open tube type of column, however, con- 
centrations vary gradually through the length of the 
apparatus and not stepwise as is the case in an ideal 
plate column. When, therefore, the data obtained re- 
quire a more fundamental interpretation, recourse 
should be had to the principles underlying the material 
transfer in the gradual countercurrent vapour-liquid 


exchange process. 
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In a recent study, Byron, Bowman and Cou t [7] 
gave a qualitative evaluation of pressure effects on the 
efficiency in the normal rectification process (“contact 
rectification’). According to their point of view, in 
the pressure region below 1-3 kN/m* (10 mm Hg) ft 
two detrimental influences of pressure are becoming 
noticeable : 

a) relative increase in vertical back-diffusion over 
transverse diffusion owing to increase in diffusion 
coefficient, 

b) non-equilibrium between liquid and vapour at 
the vapour-liquid interface owing to the relatively low 
vaporization rate. 

Both these factors would invalidate the theoretical 
approach to the vapour-liquid exchange process along 
lines customary in the case of processes at normal 
pressures. 

However, in their deductions with respect to 
point a) Byron et al. disregarded the convection factor 
in the exchange process, which counteracts the vertical 
back-diffusion. In the analysis of rectification effi- 
ciency of wetted wall columns by WestHaver [32] 
and by Kun [23] it is shown that the back-diffusion 
effect is controlled by the ratio of diffusivity to linear 
vapour velocity. Since pressure reductions affect both 
factors to about the same extent it is seen that at low 
pressure back-diffusion cannot be more detrimental 
than at high pressures. 

In discussing point b) indicated by Byron ef al. it 
should be noted that the approach to equilibrium at 
the vapour-liquid interface depends on the ratio of 
the total vaporization-condensation rate to the rate of 
material transfer. In cases where the second rate is 
large compared with the first, equilibrium cannot be 
obtained. 

Calculation according to the well-known Lane. 
MUIR formula shows that at a partial pressure of 
0-067 KN m* (0-5 mm Hg) the vaporization rate at 
100° C of a substance with molecular weight about 
200 equals about 0-2 kg m* sec. 

The transfer rate in a column as tested in the pre- 
sent investigation, though depending on relative vola- 
tility, may be estimated at about 0-0025 kg’m® sec. 
The transfer rate therefore equals about 1% of the 
vaporization-condensation rate and equilibrium be- 
tween vapour and liquid may be assumed, even at the 
low pressure quoted. 


+ In this region the mean free path of the molecules is still 
small with respect to the dimensions of the apparatus. 
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Therefore there seem to be no objections toa theo- 
retical approach based on principles which also form 
the starting point of considerations on the normal 
pressure exchange processes. 

Generally, there will be diffusional resistance to 
the transfer process both in vapour and in liquid 
phase. Using the “height of a transfer unit’ concept 
introduced by CurTon and Cotsurn [9], [11], the 
following relations may be given for the material 
transfer in the vapour phase and the liquid phase 
respectively, neglecting mutual influences of vapour 
and liquid flow: 


(H.T.U.)p 
d 


(HLT. ‘ 
. l 
= a, Rev’. 


Py 1.45 
( Ne 
1, Rev’ . Sey 


, 
S 4 
“L 


in which for the case of distillation 
Ve 


* dy 
| iZ—y* 


(H.T.U.)) 


(H.T.U.), (4) 


Generally, the liquid and vapour transfer factors 
(H.T.U.), and (H.T.U.)y cannot directly be measured, 
the interfacial concentrations x; and y; being un- 
The effect 
preted in terms of an overall vapour transfer factor, 


(H.T.Uoy, defined as 


known. overall found may be inter- 


4 
(H.T.U oy “2 


(5) 
and which is related with the separate phase (H.T.U.) 
values by 


(H.T.U)Joy = (HLT), +m) (LTUY,. 6) 


L 

The average value of the slope of the equilibrium 
line, m may be found from 

L (H.T.Uor (7 

V (H.T.U jon 7 

in which (H.T.U.)o, in analogy with (4) is given by 


z: 


(H.T.U.)o, ‘| f 


z, 


dx , 
. a ° (8) 

The coefficients a, p and q in eq. (1) and (2) being 
known, it will in principle be possible to predict the 
variations in separating power when varying the total 
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pressure on the distillation apparatus. This may be 
done by inserting the appropriate values of the phys- 
ical constants at each pressure. 

Considering the case of constant vapour and liquid 
mass flow, the influence of pressure upon (H.T.U.)p 
The 


vases are nearly constant over a W ide range of pressures 


will not be very great. ScuMipT numbers of 


and temperatures*, and therefore variations in 


(H.T.U.), originate only from variations in the vapour 
ReyNoLps number. In distillation, pressure varia- 
tions always being accompanied with temperature 
variations, the effect on (H.T.U.), found will only be 
due to the temperature influence on the Vapour \ Iscos- 


itv. This influence is only small (4, being approxim- 


ately proportional to | T). In accordance with these 


considerations, GILLILAND and SHEeRWwoop | 18] found 
in vaporization experiments in a wetted wall column 
that under (abt 
IS kN m*—260 KN m®*) (H.T.U.))» remained essentially 


constant at constant vapour mass flow. 


wide variation of the pressure 


On the other hand, the temperature variations 
cause large variations in the liquid viscosity and dif 


fusivitv. Using the relation 


D, My 


constant, 
1 ‘ s 


it is found that the liquid ScumiprT number is pro- 


> + 


portional to «) To,. while the liquid ReyNo ips 


number varies proportionally with «,'. Though de 
pending numerical values of the 


upon the power 


coefficients p, and q, in eq. (2), 


in general the result 
will be that a reduction in pressure causes an increase 
in the (H.T.U.), value. 

Therefore, if the liquid phase resistance cannot be 
neglected, it is seen from eq. (6) that the overall 
separating power of a distillation column will decrease 
with decreasing total pressure (m and L | remaining 
constant). 

In order to arrive at a more quantitative interpreta- 
tion of the effects described above the coefficients a. 
pand q in eq. (1) and (2) should be evaluated. Though 
this has to be done at least partly by experiment 
some general remarks can be made. 

When the vapour flow is laminar, the distilling 
column being a eyvlindrical tube, eq. (1) becomes [32]. 
(H. > mn" WODTS Rey Sep. (9) 


* These considerations are only valid when the pressure 
in the distillation apparatus ts high enough for the mean free 
path of the molecules to be small with respect to the dimensions 
of the apparatus. In the normal countercurrent apparatus 


this will always be the 


Cast. 
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Under turbulent flow conditions in the vapour, as 
will generally be the case in a Vigreux column, the 
exponential values py and q, will differ largely from 
unity. Recent investigations into the influence of 
REYNOLDS number and ScumipT number over wide 
ranges upon material transfer in open tubes and 
packed beds suggest that p, will be equal to about 
(2 [22], [27] and that gp may be put at 0-6-0-8 [21], 
[22], [29]. 

For the liquid phase when flowing in laminar 
motion along the wall of the distilling column. the 


following relation holds [35]: 


(H.T.U.),, 


a 4 , , 
_ O-O5S89 Re, Se, . (10) 


Practical application of eq. (10) is complicated by the 
fact that the liquid layer thickness, w, depends on 
the liquid flow rate. In an open tube column some 
indication about the magnitude of this variable may 
be found from hold up measurements. 

Generally, in distillation, full turbulent flow does 
not develop in the liquid phase: very little is known 
about the values to be assigned to the coefficients p; 
of some values of 


and gq, in this case. Recalculation 


BrinsMaADE and Butss [6], who measured transfer 


coefficients in the extraction of acetic acid from 
methyl isobutvl ketone with water in a wetted wall 
column, shows a value of about 0-65 for q,: the value 
of p, equalled zero, indicating that the liquid flow 
rate had no influence. 

From the above considerations it may be concluded 
that unless the liquid and vapour phases are flowing 
in Viscous motion, most of the coefficients in eq. (1) 
and (2) have to be assessed by experiment. An ex- 
ception may be made for the power coefficient qy 
on the vapour ScumMipT number. The actual variation 
number being small, the value 


in this ScHMIpDT 


to be used for gy does not matter very much. 
In this work gy, will be put at 0-67 without further 
verification. 

The evaluation of the remaining constants should 
be done by calculating (H.T.U.), and (H.T.U), values 
from the overall (H.T.U.)y, values. This can be done 
with the aid of eq. (6) and (7) when data from ex- 
periments in which m or LV have been varied are 
available. 

Some final remarks should be made about the in- 
fluence of mechanical agitation, as is applied in a rotat- 
ing band column. By this agitation material transfer 


is increased by creating extra turbulence in the vapour 





F. J. Zctperwee: Vacuum distillation I 


phase particularly. As a first approximation the in- 
fluence of Vapour velocity and extra turbulence in- 
duced by the rotation on (H.T.U.),; may be visualized 
as separate independent effects*. Generally, this 
will show up in a decrease of the coefficient a) in 


eq. (1) as a function of the intensity of agitation 


From the studies of TayLor [30] on the stability 
of a stagnant fluid between two rotating cvlinders it is 
shown that turbulence only 
sets in when the rotational 
speed exceeds a certain 
critical value. The existence 
of a critical rotational 
speed in the case of a fluid 
flowing in axial direction 
between the evlinders was 
shown by Cornisn [12 
by pressure drop measure- 
ments on the water flow 
through narrow annuli with 
rotation of the inner ey 


linder 


It should therefore be 
expected that in the case 
of rotation of a strip ora 
eviindrical tube there will 
also be a critical rotational 
speed below which rotation 
does not increase the ma- 
terial transfer effectively. 
In a first approximation 


assuming the width of the 





Fig. 1. Photograph of part strip to be equivalent to 


of Vigreux column tested the diameter of a rotating 
evlinder, a relation for the 
critical ReyNoips number of rotation may be de 

rived from the equations given by Face [14 

ordea d 1 
(Re »,) zt ~ (11) 
. My 1 

It will be clear that eq. (11) cannot be exact for 
a rotating band column: a distinct influence of the 
dimensions of the apparatus on the critical rotational 


speed is to be expected however 


* In this approximation it is assumed that the resi 
dence time of the Vapour entering the bottom of the rotary 
column is long enough to ensure the vapour obtamung the 
flow puittern induced by the rotating element This will only 
be the case when a relatively long (° mfimuite column is 


bemg used 
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The above discussions suggest that the reduction 
in coefficient ay in eq. (1) under influence of rotation 
should be correlated against the ReyNoLps number 
of rotation, Re;,. As shown by studies on heat and 
mass transfer in agitated vessels [10], [19], [20]. the 
power coefficients of the Scumipr number may be 
put at 5 7 when turbulence has been fully develop 
ed: as before, a value of a (67 will be assumed 
in the present study. 

Actually, the effect of agitation will only be great 
in cases where vapour flow is so small that laminar 
motion prevails when the stirring device is not ac 
tuated. In this case (H.T.U.), will increase linearly 
with Re, 

The ultimate correlation equation by which the 
vapour phase transfer factor should be represented 
will therefore have the following general appearance 

(H.T.U.)) 


d ; Re ti Re i Ncj a 


a R 


APPARATUS AND EXPERIMENTAL PROCEDURE 
Viare ux column 
The measurements with the Vigreux column may be 
divided into two parts, ¢.9 measurements of separating 
power, pressure drop and maximum boil-up rates at 
pressures from about IS KN m*® 101-5 KN m*® (1 to 
760mm Hg) on the one hand and determination 
of the separate (H.T.U.), and (H.T.U.), values at 
atmospheric pressure on the other hand 

The low pressure performance data were obtained 
in a Vigreux column with a rectifying section of 1m 
in length and 11 mm in diameter. The vertical distance 
between two rows of horizontal indentions was 24 mm 
A photograph of a part of the rectification section 
used is given in Fig.l. The column was surrounded 
by a heating jacket to secure adiabatic operation, 
Boil-up rate was measured at the bottom end of the 
column by means of a calibrated measuring tube in 
the reflux return line. To prevent irregular boiling 
and bumping, the contents of the boiling flask were 
heated by means of a bare nichromium wire coil. To 
measure the pressure drop an oil differential mano. 
meter was fitted between the boiling flask and the top 
of the condenser of the column. Some separate press- 
ure drop measurements were made by recording the 
temperature drop between top and bottom of the 
column when distilling a pure liquid (¢.g. n-hexadecane ) 
at total reflux. 

The separating power of the column described was 


determined at atmospheric pressure (”-heptane methyl- 
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cyclohexane, n-decane, transdecalin, Tables 1 and 2), 
at 13-3 kN m? (100 mm Hg) (n-decane transdecalin, 


Table 1. 
Binary system n-heptane-meth yleyclohe cane 
101-5 kN /mé 


Separating power of Im Vigreux column 


Pressure 





Boil-up 
rate 


nae nee Theoret- 
D D ical plate a 


ml h top bottom n 
1-4098 
1-4089 
1-4088 
1-4088 
1-4087 
1-4087 
14005 
1-4094 
1-4090 
1-4088 
1-4091 


-4030 
4009 
-4020 
4007 
4018 
“4015 
4021 
4026 
4018 
4028 
4011 


740 
410 
0 
220) 
520 
+0) 
240) 
Ho 
610 
S70 
220 


1050) (max) 





Table 2. 


Binary system n-decane-transdecalin 


Separating power of lm Vigreux column 





Theoretical 
plates 


Mol. 


% n-decane 


Boil-up rate 


Pressure 


kN me ml h top hottom 


1O1-5 930) 
L050 (max) 
340 
280) 
7) 


70-3 


65-1 


380) 
250 
140 
240 
450 
370 
SOO (max) 
19) 
ow 
250 
160 
280) 
230) 
190 


350 (max) 





Table 2), at 1-33 kN /m? (10 mm Hg) (n-decane, trans- 
decalin, n-hexadecane/n-heptylbenzoate, Tables 2 and 


14° 
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3), and at 0-16kN/m* (1-2 mm Hg) (n-hexadecane! 
n-heptylbenzoate, Table 3). 

Samples were taken at the top and bottom of the 
column, using special vacuum adapters. An equilibra- 
tion time of 30 to 60 min was allowed for. Analysis 
of the products was made by determination of re- 
fractive indices with a normal ABBE refractometer, 
read off to the fourth decimal place. The n-heptane- 
methyleyclohexane data were converted directly into 


Table 3. 


Binary system n-he xadecane-n-he pt ylbenzoate 


Separating power of Im Vigreux column 





Theoretical 
plates 


Pressure 


Boil-up 


Mol. % nt 416 
rate drop 


Pressure 


LN om? mlh kN im? top bottom 
1-33 170 
20) 
40) 
410 


380) 


0-09 74-4 1-9 
74-4 1-9 
0-26 76-0 4-6 
O30 76-0 4-6 
0-30 75-7 4-6 
170 0-09 


21 


450 (max) 

270 45 
200) 0-30 
230 0-39 
170 0-26 
120 O17 
120 O13 


280 (max) 





plate numbers with the aid of the alignment chart 
given by Brooks, NELSEN and Zaun [5]; refractive 
index data from the other binary mixtures were 
converted into plate numbers using the calibration 
curves and relative volatility data given in the com- 
panion paper [36]. Maximum boil-up rate was record. 
ed when initial flooding was observed visually at some 
point along the length of the column. 

As described in the theoretical part, a more tho- 
rough interpretation of the results to be obtained at 
various pressures may only be given in terms of the 
separate phase resistances. This may be done by 
analyzing the overall performance data, expressed as 
(H.T.U.joy, by means of eq. (6), using data obtained 
at various m and L/V values. In order to reach some 
accuracy in these evaluations, the variations in m or 


LV should be as large as possible. This implies that, 


in the case of distillation, a binary mixture should be 
used whose components have widely different boiling 
points. Further, variation in m can only be obtained 
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when the average composition of the contents of the 
column is varied: this can be done if the separating 
power of the apparatus is limited. 


Some further points should be considered when 
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n-octane system, the concentrations in the lower part 
might be expected to be constant and to be equal to 
This was checked by tak- 


ing samples at points at about 0-10 and 0-20 m from 


the “pinch composition ”’. 





























evaluating the separate phase resistances by this the bottom end of the packed section. 
method. The results are reliable only when + 
under the conditions applied, no channelling = 
of the liquid phase or incomplete wetting of g “i fh 
the column area take place. These effects may a ]J 
give rise to discrepancies, which cannot be ae 
correlated by eq. (6), as may for instance be zz 
shown for data of Duncan, Korrout and Wrrn $ _ rad] { COTTON WOOL 
row [13] obtained on a Rascuic ring packed + | 
column. Moreover, when varying the LV ratio 
the possibilities of changes in the vapour-liquid ali ; = VIGREUX SECTION 
interaction should not be overlooked. XN 

In the present work it was found that only ro Zid >= TOP PRODUCT LINE 
very few wide boiling binary systems do not JZ q ¢ 2mm 
exhibit the channelling effect when distilled in HERMOCOUPLES 
an open tube type of column. For this reason ae i 
use was made of the binary mixture of 2.2.4 HEATING JACKET. r akon oleae 
trimethyvlpentane-n-octane 

The experiments in which m was varied L 
(at total reflux) were made in a short (0-22 m) j a 
Vigreux column, cut from the l-m column de- 
scribed before. This column was again provided [ 
with reboiler, condenser and reflux-measuring BOTTOM “ PINCH’ SAMPLES 
device; adiabatic conditions were maintained SAMPLE 
with a heating jacket, in which the heat ap REFLUX - 

, a . MEASURING 

plied was controlled by a lagliabue Celectray TRAP ,7OP SAMPLE, TOP PRODUCT 
which kept the temperature difference between PUMP 
column surface and heating jacket at zero eae 
value. REBOILER 

Variations in LV ratio were obtained by ' 
making experiments at partial reflux. In order 
to determine the exact liquid-vapour ratio 
use was made of an experimental procedure 
in which two points of the operating lime were Fig. 2. Dinguam of eppesatus fer niaueniia a ii eile 


determined, viz. the distillate composition and 

the liquid composition in the intersection point 

of equilibrium line and operating line (“pinch com- 
position"). This was achieved with the experimen- 
» 


tal lay out, sketched in Fig. 2. 
The small length (0-22 m) of Vigreux column was 


of 


a packed column, filled with small nichromium wire 


placed upon a 0-50 m section, diameter 16 mm 


rings (about 2. mm in diameter). The separating power 
of this piece of packed column was equivalent to 


about 20 theoretical plates, so that under conditions 
@@ 


of partial reflux, using the 2,2,4-trimethylpentane- 


ratio with Vigreux column 


The top end of the Vigreux section was fused to 


a short wetted-wall column (020m) in order to 
ensure that the reflux entering the Vigreux section 
was at boiling point. The top product was with- 
drawn at the top of the Vigreux section by means 
of a small reciprocating pump with variable stroke. 
The distillate heated to boiling point, was delivered 
continuously to the reboiler, thus ensuring a closed 
system in which the concentrations were becoming 


stationary. 
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Adiabatic operation of the apparatus was again 
effected by means of automatic control of the heating 
jacket with the Tagliabue Celectray; two control 
points and a heating jacket with two separate windings 
were now used, one at the Vigreux section and one at 
the packed “pinch” section. 

Apart from the samples withdrawn from the 
pinch section mentioned above, liquid samples were 
taken at the bottom end of the Vigreux section and of 
the distillate withdrawn with the pump. All samples 
were withdrawn continuously by means of small 
capillaries at rates of 1-2 ml! per hour. In each ex- 
periment the rate of distillate was determined separa- 
tely, from which by means of the proper LV ratio 
(from top and * pinch” compositions) boil-up rate and 
reflux rate were calculated. Generally there was good 
correlation between the calculated reflux rate and the 
reflux rate measured at the bottom end of the packed 
section with a reflux measuring trap. 

Analysis of the liquid samples was made by deter- 
mination of the refractive index with a Bausch and 
Lomb precision refractometer, enabling readings in 
three units in the fifth decimal place. By this means 
concentrations could be determined with an accuracy 
of about 0-5 mol % absolute. To ensure utmost pre- 
cision, the seale of the refractometer was calibrated 
with pure n-octane or 2,2.4-trimethylpentane during 
each determination. 

From the concentration data obtained, the separat- 
ing power of the Vigreux column was calculated, as- 
suming the 2,2,.4-trimethylpentane-n-octane system to 
be ideal; corrections were made for deviations from 
ideal gas laws [1]. The calculated relative volatility 


data, based upon vapour pressure data from A.P.I. 


project 44, National Bureau of Standards, Washington, 
are given in Table 4. Calculation of the total reflux 
data was made by using a differential liquid and 


The 
partial reflux experiments were calculated by plotting 


vapour transfer units vs composition curve. 


the top, bottom and pinch concentrations on a large- 
scale McCape-THIELE diagram, drawing the operating 


Table 4. 


2,2,4-trimeth yl pentane-n-octane 


Calculated relative volatility 
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line through the top and pinch compositions and read. 
ing (y* —y) and (x 
respectively. (H.T.U.)oy and (H.T.U.)9, were found 
by graphical integration of the results obtained. The 


x*) values as functions of yand x. 


average slope of the equilibrium line was calculated 


from eq. (7), using LV values read off from the 


McCape-THIELE diagrams. 


Rotating band columns 


The rotating band columns were purchased from 
H.S. Martry, Evanston, Ill., U.S.A. The rectifying 
section of the first column, (a ** Pyros-Glover micro- 
still”), in which the majority of the experiments were 
run, consisted of a straight glass tube, 0-60 m in length, 
6mm internal diameter. In this tube a flattened coil 
of nichrome wire (abt. 0-2 mm thick) through which 
a nichrome strip is inserted could rotate. The width 
of the strip was about 4-5 mm, the total width includ- 
ing the nichrome wire about 5mm. The strip was 
rotated by means of a variable speed motor, driven 
via a tungsten shaft, running through the condenser 
section. The second * Pyros-Glover microstill” tested 
had a rectifying section of 3 mm diameter at a length 
of 0-40 m. The total width of the rotating band in this 
column was about 2mm. The microstills were fitted 
with a device for partial condensation (the primary 
condenser cooled with air, the secondary or product 
further 


completely surrounded by a vacuum jacket and heat- 


condenser with water). The columns were 


ing jacket. In order to avoid bumping in the low 
pressure determinations, heating of the reboiler con- 
tents was again done by means of a bare heating 
spiral. To minimize effects of draft upon the boil-up 
rate, reboiler and lower end of the columns were placed 
in a Dewar flask. 

In the low-pressure experiments the inlet of the 
shaft of the rotating strip should be sealed; this was 
reached by piercing the shaft through a cylindrical 
piece of rubber, which was fitted to the top of the 
condenser by means of a thick vacuum tube, glycerol 
being used as a lubricant. 

In determining the separating power of the micro- 
stills, samples of the top product were drawn by re- 
ducing the cooling of the primary condenser somewhat 
after an equilibration time of 1-2 hrs. No samples of 
the bottom contents were taken during the distillation : 
the bottom composition was analyzed before and at 
the end of a distillation and the average value taken. 
The hold up of the column being small, this method 
was believed not to introduce serious errors. 
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The boil-up rate was measured by counting drops 
from a calibrated dripping tip at the bottom end of the 
rectification section, 

Experiments were run at atmospheric pressure in 
the 6mm and 3mm stills with »-heptane-methyl- 
cyclohexane, varying boil-up rate and speed 
of rotation (Table 10). Reduced pressure ex- 
periments were made in the 6mm still with 
10-7 KN m® and at 


with 


n-decane-transdecalin at 
24) KN m® (80 
n-hexadecane-n-heptylbenzoate at 0-67 KN m? 


and 15mm Hg) and 
(5mm Hg): the rotational speed was kept 


constant at 50 rev. p. second (Table 11). 


Discussion oF EXPERIMENTAL RESULTS 
Vigreux column 
The results obtained with the l-m Vigreux 


column at various pressures (Tables 1, 2, 3) 


— ——. —4 aCz7- mech. 
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nearly constant at 9-10 plates. At constant boil-up 
rate, however, a distinct pressure influence is shown: 
decreasing the pressure from atmospheric to 1-3kN/m? 
(10 mm Hg) shows a decrease in separating efficiency 


25%. 


of about This effect, as already described in 


THEORETICAL PLATES 


‘LEGEND + T T | ] 
—4——-+ MCyQ-ft.dec. —+—+MCqg-nhb 


, + KN/m2?_, i +— KN/m2_| i 1 KN/m2 


1015 101.5 
t t 1993-7 13 
13 » 0.16 


i 








are plotted as theoretical plate number rs 
boil-up rate in Fig. 3. It is clear that the 
separating power of the column is not very 
much influenced by pressure. Comparing the plate 


number at maximum boil-up rate, efficiency is 


Table 5. Separating power of 0-22 m Vigreux column 
(total reflux) 
oo 


Binary system 2,2,4-trimeth yl pentane-n-octan 


101-5 kN ni 


Pressure 





Boil-up rate ( ‘om position, mol. fr. Separating 


power 

ml h top hottom H.T.U Wo, 
141 
OF 106 
OLS 


loa 47-0) 
420 a2 
630 MS 
loon 46-0 
S40 45-5 
420 52-0 
690 47- 

1a 41- 135 
lho 112 
7 is . 005 
420 f 33- 0-087 
blow O-113 
420 s-5 : 0-OSS 


1000 0-097 
1100 0-108 
420 s2- . 0-080 
1050 
420 “ Sle 0-080 
420 . . 0-082 
720 70: 0-091 





200 400 600 800 1000 
BOIL - UP RATE, mi/br 


Fig. 3. Separating power of Lm Vigreux column at various pressures. 


the theoretical part, may perhaps be due to an in- 
crease in liquid phase resistance at lower pressures, 
In order to arrive at a more quantitative inter- 
pretation, separate phase resistance determinations 


(nT), (m) 
0 16 


0% = = Rey 3100 
Rey 2200 


Re, 1100 





HT 
wih My 





0.042 
0.047 
0.052 




















a2 04 06 08 90 92 J 
im. SLOPE OF EQUILIBRIUM LINE 








Fig. 4. Determination separate phase resistances of 


Vigreux column. 


were made at atmospheric pressure with the 0-22-m 
Vigreux column already described. 

Results at total reflux, in which the average con- 
centration in the column and thus m was varied, are 
listed in Table 5. From these data a plot of (H.T.U.)oy 
vs boil-up rate was made; at three distinct values of 
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the boil-up rate (eg. at 400, 800 and 1100 mlhr) 
(H.T.U.)oy data were read off and plotted in Fig. 4 
as a function of the average slope of equilibrium line, m. 
In this plot, the intersection points on the ordinate 
the (H.T.U.), values, 
the slope of the lines representing the (H.T.U.), 


indicate—according to eq. (6) 
values. From the data obtained (recorded in Fig. 4 
too) it is seen that the vapour and liquid phase resis- 
tances are of the same order of magnitude, liquid 
phase resistance being the largest of the two. 
TU), (m) 


, 


+ TOTAL REFLUX (VARYING VAPOUR FLOW) 





— 





20 40 60 80 100 
RE YNOLOS NUMBER OF LIQUID, Re, 
Liquid phase resistance of Vigreux column at 


total reflux and partial reflux. 


Fig. 5. 


The (H.T.U.)y data were used to evaluate the 
coefficients ay and py still unknown in eq. (1). 

The following equation was obtained *, the diam. 
eter of the column being included in the constant ay: 


(H.T.U.), 


Calculated and experimental (H.T.U.), data are 


= 0.0146 Ref Sci” . (12) 


compared as follows: 





(A.T Uy, (m) (ALT Uy, (m) 


Rey eq. (12) 


erp. 


1100 
2200 
3100 


0-042 
0-047 
0-052 


0-0421 
0-0482 
0-0519 





* The physical constants at distillation conditions necess- 
ary for calculation of the Re and Se numbers for the 2,2,4- 
trimethylpentane-n-octane systern and the other systems are 
given in Table 8. Liquid densities and vapour viscosities were 
taken from current handbooks like 
Tables ete.; vapour densities were calculated from ideal gas 
law, liquid viscosities were estimated from the work of NepER- 
BRAGT and BorELHOUWER [24]. Finally, liquid diffusivities 
and vapour diffusivities were calculated by the correlations 
given by Wiike [33] and Griuimanp [17], respectively. 


International Critical 
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be concluded that full turbulence has 


It 


developed at REYNOLDS numbers as low as 1100; this 


may 


is caused by the disturbances induced in the vapour 
by the indentions in the column surface. 

The liquid phase transfer factors (H.T.U.), are 
plotted in Fig. 5 as a function of ReyNoLDs number 
of the liquid. It is seen that (H.T.U.), increases with 
Re,,, raised to a power of about 0-3. This indicates 
the 
REYNOLDs numbers are very low, nearly all below 
In that 
case, (H.T.U.),; should increase about linearly with 
Re L> 


pected that in the measurements under total reflux 


turbulent flow in the liquid phase. However, 


100, and laminar flow should be expected. 


as indicated by eq. (10). It was therefore ex- 
turbulence was induced in the liquid phase by the 
vapour flow. Indications of this effect were already 
the 
l-m column. At boil-up rates much lower than the 


found when making the low pressure runs in 
maximum boil-up rate the liquid phase was visibly 
agitated by the vapour flow. At higher vapour velo- 
cities the turbulence induced in the liquid phase will 
This effects a simultaneous decrease in the 
(H.T.U.), does not in- 


crease so sharply as would be expected on the basis 


increase. 


liquid phase resistance, /.¢. 


of pure laminar behaviour only. 

In order to evaluate this effect further, experi- 
ments at partial reflux at nearly constant vapour rate 
were made. In these experiments, the turbulence 
effect of the vapour flow on the liquid phase may be 
assumed to be nearly constant. The data obtained by 
the method described in the experimental part are 
assembled in Table 6. In order to evaluate the sep- 
arate (H.T.U.),; values from the (H.T.U.)oy values 
measured, use was made again of eq. (6), in which the 
vapour phase transfer factor (H.T.U.), was calculated 
with eq. (12). The (H.T.U.), values ultimately obtain- 
ed are plotted against ReypoLps number in Fig. 5 
too. 

Though the resulting data are of only moderate 
accuracy, it is clear that at a constant vapour velocity, 
(H.T.U.), increases about linearly with Re, , i.¢. the 
separating efficiency is proportional to contact time 
as is the case in pure laminar flow [eq. (10)]. From 
Fig. 5 and Table 6 it can be deduced that at constant 
liquid rate doubling of the vapour rate effects a con- 
siderable decrease of the liquid phase resistance 
[compare the (H.T.U.), values under total and partial 
reflux at Re, about 60}. 

The results thus far obtained on the short Vigreux 
column have given the following information: 1) a 


} 
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Table 6. Separating power of 0-22 m. Vigreux column at partial reflux 


Binary system 2 


,2,4-trimeth ylpentane-n-octane ; pressure 101-5 kN m* 





Boil-up 
rate 


Liquid 
rate 


ml h ml h top hottom 


1380 720 0-220 
O-314 
(455 
0-470 


0-185 


450 
680 0-623 
620 


1o50 
1280 
1280 
1310 
1180 


0-730 
0-425 
475 


0-556 


970 
60 
TO 
1250 


0-205 
1220 (1-265 


1250 


Compositions mol. fract. 


(H.T.U Joy (H.T.U) yp (H.7T.U.), 


pinch m m 


195 
(230 
0-400 
0-405 
OLLO 


(144 
0-110 
0-118 
0-123 
0-139 
0-129 


0-054 
O-O51 
0-053 


0-037 
0-036 
0-036 
0-040 
0-048 
0-049 
0-031 
0-066 


0-053 
0-053 
0-095 0-052 


215 O-1LLO 0-052 


0-053 





relation for calculation of vapour phase transfer factors 
(H.T.U.)y [eq. (12)]; 2) liquid phase resistance is of 
the same order of magnitude as vapour phase resis- 
tance at atmospheric pressure ; 3) liquid phase transfer 
resistance is decreased by turbulence, induced in the 
liquid by the vapour flow. 

This information should be used for interpretation 
of the data found at various pressures with the l-m 
the effect 
under 3) greatly complicates such a further inter- 


Vigreux column. However, mentioned 
pretation. No adequate data at the various pressures 
being available at present, the following approximate 
reasoning as to the turbulence effects in the liquid 
phase due to vapour flow may be employed. 

The action of the vapour flow will be a result of the 
impact and shearing forces exerted by the vapour on 
the liquid phase. The degree to which the liquid phase 
is being agitated by these forces depends not only 
upon their magnitudes but also on the liquid flow rate 
and physical properties of the liquid. These same fac- 
tors however play a similar role in the flooding pheno- 
menon in which the liquid layer is disrupted and 
thrown upwards by the forces exerted by the vapour 
flow. Therefore, as a first approximation, it seems 
plausible to take the turbulence effects in the liquid 
phase in the total reflux experiments the same at 
equal fractional flooding rates. Thus, column efficien- 
cies at various pressures should be compared for in- 
stance at maximum boil-up rate or at half the maxi- 
mum boil-up rate. Under these conditions of constant 
turbulence effects, (H.T.U.), will be assumed to 
increase linearly with Re, as found with the octane ex- 
periments at atmospheric pressure (Fig. 5). 

The above considerations enable the calculation 
of (H.T.U.), data from the overall separating power 


figures obtained at various pressures in the l-m Vigreux 
column. In using eq. (6), (H.T.U.)p» values were first 
calculated with relation (12) from the octanes experi- 
ments. According to eq. (2), the results should ulti- 
mately be correlated against the liquid Scumipr 
number, using (H.T.U.), Re; as the proper variable. 

Calculations of (H.T.U.), Re, values at maximum 
boil-up rates and half maximum boil-up rates are 
summarized in Table 7. The results are finally plotted 
in Fig. 6 vs the liquid Scumipr number. Correlation 
of the data by this method is fairly good; as might be 
expected two separate functions are obtained, namely: 


At maximum boil-up rate 


(H.T.U.), = 12-4 « 10~* Re; Sef’. (13) 


At half the maximum boil-up rate 


(H.T.U.), 


11-8 « 10 ® Re, Ne,. (14) 


Fig. 6 shows that at maximum boil-up rate 
(H.T.U.), Re, values are about half the values obtain- 
ed at half the maximum boil-up rates. This indicates 
the effect of higher turbulence in the liquid phase at 
the higher vapour rates, which is reflected in power 
coefficient of the Scumipt number in eq. (13), which 
differs from unity. 

At half the maximum load however this power 
coefficient equals unity; the turbulence effects of 
vapour flow are obviously negligible now. 

An explanation may be given by the consideration 
that at half maximum load the shearing forces exerted 
on the liquid are only about a fourth of the stresses 
at maximum boil-up rate (taking the friction factor 
approximately constant). 

Assuming eq. (14) as a true expression for the 
laminar flow condition, the ratio of the constants in 


184 
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Table 7. Calculation of (H.T.U.);, in Im Vigreux column at various pressures 





Boil-up 
rate 


System ml h 


half 


kN om? mar. 


mar. 
OLLO 3200 
1600 
3000 


lO}. 0-109 
0-091 


0-102 


n-Cy-meh 1050 
0-092 
1OL-5 0-107 


nr t-clec 1O50 


0-087 0-092 1500 


0-108 0-109 
0-093 
O-114 


O-105 


13+: SOO 2100 
1050 


0-008 
120 
O-lil 
0-099 


850 


450) O-1L07 2200 


0-098 0-090 1100 


O16 2so0 0-109 O- 100 


0-108 O-005 
2,2,4-i-m-p- LOLS 


nC, 


L500tt 
THOtt 


yh 


m(H.T.U. 


0-051 . . S: lo 
0-044 . . 1o~ 
0-046 lo 
0-040 . . 1oO7 
0-044 . , . 22 ?- 10 
0-038 : . . “{ lo” 
0-040 . Se OF { 7-4 lo-: 
0-035 , lt . - 5 lo” 
0-040 . “Ss ‘OF 3°: 1o7 
0-035 “O55 Ss: . 6 lon 
0-037 . S: . 10 
0-032 . “S! . 7 lo~ 
10 





* Data taken from Fig. 3. + Data from Fig. 5. 


** Calculated with eq. ( 12). 


eq. (14) and in eq. (10) should give an indication of 
the average liquid layer thickness. This ratio shows 


(mTu) 
—_ 


AT HALF 


s LEGEND 
FLOODING RATE,” 
/ 


SYSTEM 








octanes 
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Fig. 6. Correlation of liquid phase resistance data of 1 m 
Vigreux column. 


for w a value of 0-20 mm. Though no special hold up 


measurements were made with the present column, its 


tt Caleulated on average column pressure and temperature 


value may be estimated at about 6-7 ml [8]; this 
shows an average liquid layer thickness of 0-17 to 
020mm. The agreement between the two w-values 
thus calculated is quite good and in fact much better 
than would be expected on the basis of the many 
assumptions and uncertainties in physical constants 
involved. 

It has been shown that the method of interpretation 
used leads to quite logical and consistent results. 
They enable us to formulate the following insight 
gained with respect to the influence of pressure upon 
the separative effect of a Vigreux column. 

At constant boil-up rate, a decrease in pressure 
gives rise to a large increase in liquid phase resistance 
to material transfer because of increasing viscosity. 
As at atmospheric pressure this liquid phase resistance 
is already more than half the total resistance, it causes 
the overall separating power to diminish at a reduction 
in pressure. However, this effect is partly nullified by 
a simultaneous increase in turbulence induced in the 
liquid by the increasing vapour velocity. 

The question may be raised how far the present 
results may be compared with those obtained on 
other types of column. 

In the case of packed columns, the effect of vapour 
rate in generating turbulence in the liquid may be 
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much smaller than with the Vigreux column. This 
will be caused by the distribution of the liquid over 
more or less capillary spaces in which it cannot well 
be agitated, in particular when the packing material 
is fine. Moreover, wetting of the packing may be 
a function of the liquid rate, complicating the (H.T.U.), 
vs flow rate relation. Both factors suggest that in 
packed columns a decrease in pressure will cause 
a greater loss in overall efficiency than in the Vigreux 
column. 

This reasoning is substantiated by available lit- 
erature data. Fretpman ef al. [16] showed a Pod. 
bielniak Heligrid column to lose about 50° of its 
separating power when reducing the pressure from 
101-3 KN m* to 6-7 KN m* (760 to 50 mm Hg). Approx. 
imately the same effect may be shown from data 
obtained by Struck and Kinney [28] on several 
types of column packing. It should be kept in mind 
however that in their work incorrect values of the 
relative volatility for the n-decane-transdecalin system 
were used. Correction of the Srruck and Kinney 
data using the relative volatilities mentioned in the 
present paper shows a decrease of about 30% in sep- 
arating power at the same pressure region mentioned 
above. 

Bere and Porovac [2] determined the efficiency 
of Fenske helices packing at various pressures with the 
toluene-n-octane system. The vapour-liquid equili- 
brium data used are thermodynamically inconsistent 
however. Approximate correction of these data shows 
that at maximum boil-up rate the separating efficiency 


decreases by 20% when applying a pressure reduction 


from atmospheric pressure to 6-7 KN m* (50 mm Hg). 


Table 8. 
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Apart from the study of the pressure effects upon 
the separating power, some other aspects of a frac- 
tionating column should be considered in judging the 
performance of the apparatus at low pressures. 

Pressure drop data obtained with the n-hexadecane- 
n-heptylbenzoate system are given in Table 3. It is 
seen that at the lowest pressure used (0-16 kN m? or 
1-3 mm Hg) pressure drop of the Vigreux column is 
still moderate though equal to a number of times the 
top pressure (at max. boil-up rate about 0-4 kN m?* 
(3mm Hg). Reep and Fenske [25] suggest that 
pressure drop data of a packed column obtained at 
various pressures may be correlated by plotting the 
product of average vapour density and pressure drop 
against mass flow rate. As demonstrated in Fig. 7, 
this correlation also holds fairly well for the Vigreux 
column. 

The pressure drop measurements recorded with 
the oil manometer being taken between reboiler and 
top of the column, it was thought that some extra 
pressure loss occurred in the reboiler section (including 
the reflux measuring device). Therefore some pressure 
drop measurements on the Vigreux section proper were 
made by recording the temperature drop between top 
and bottom ends of this section when refluxing pure 
n-hexadecane (see Table 9). The temperature drop 
data were converted into pressure drop data by using 
the vapour pressure curve of n-hexadecane. The re- 
sults (plotted in Fig. 7 too) show that the actual 
pressure drop over the Vigreux section is some 30% 
lower than the pressure drop over the whole column. 

Maximum boil-up rate is seen in Fig. 3 to be af- 
fected markedly by a reduction in total pressure. As 


Constants used in calculations 





Liquid 


Pre assure 


S yatem De naily Viscosity 


kN m* kg ni N sec/m? 
n-heptane m-cyclo- 
hexane 
2.2,4-1-m-pentane- 
n-octane . 
n-decane-transdecalin 
n-decane-transdecalin 
n-decane-transdecalin 
n-hexadecane- 
n-heptylbenzoate . 
n-hexadecane- 


n-hept ylbenzoate 


Difhusivity 


m?* sec 


Vapour 


Density Viscosity Diffusivity 


kg in N sec m?* m?® sec 
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a first approximation it may be concluded that by approximate statement is in agreement with the results 
employing a tenfold reduction in pressure, maximum found by Struck and Krxney [28] on packed columns. 
boil-up rate is reduced to half its original value. This A more accurate correlation may be given by plotting 
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Fig. 8. Correlation maximum throughput data of Vigreux 
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rate) versus the ratio of vapour and liquid den- 
sities [25], as shown in Fig. 8. 








' 
N 

10 16° 1 - “am — 
FLOW RATE, kg/sec 50 a  TEGENO 





Fig. 7. Correlation of pressure drop data obtained with | _|SY¥MBOL| SYSTEM 
the n-hexadecane n-heptylbenzoate system in Vigreux | ' 
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Table 9. Pressure drop measurements on Vigreux column, 

















length 1-0 m, diameter 11 mm, at 0-16 N/m* top pressure by 


means of temperature drop when refluxing n-hexadecane 





Vapour 


Boil-up Pressure 


temperatures 
rate drop 
° ( ’ 


ml/h top bottom kN | 


270 109-5 130-5 0-40 
220 110-0 =: 128-0 0-33 
190 109-0 =: 124-0 0-24 
150 109-0 —s-«121-0 0-17 
120 109-0 118-0 0-12 _—-_ -_ = x: m 

160 109-0 120-0 O15 mi/hr 

200 109-0 123-0 0-22 Fig. 9. Separating power of 6 mm spinning band still at 
260 109-0 125-0 0-26 various pressures. 

300 109-0 129-0 0-36 

310 110-0 = 132-0 0-44 No measurements of hold up were made in this 
150 110-0 «1225 | = 0-20 study of the Vigreux column. It is believed that its 


120 109-0 115-0 0-07 value will not be much influenced by the pressure 
240 108-0 =—-:127-0 0-32 
250 109-0 133-0 0-49 
140 108-5 118-0 0-13 








applied, viscosity and density effects entering the 
equations for layer thickness of a falling liquid film 
only by the 1/3 power [15]. 
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From the results obtained in this Table 10. Performance of rotating band columns at atmospheric pressure 


study it may be stated in conclu- System n-he plane-m-c yclohexane 





sion that, though of only moderate Rotational 


Column Boil-up rate nD products Separating power 


separating power, the Vigreux type yee speed 
of fractionating column will be a (H.T.U Joy 
. . . . : : 7 , " J ‘ , ? he “” 

useful tool in vacuum distillation _ mp | Boy |sccjers) Bop ~~ onus 


down to bottom pressures of about 
0-67 KN m? (5 mm Hg). 1950 
2000 


3963 
3975 
4000 
007 
“S095 
4055 
4075 
4102 
4105 
4117 
“3987 4165 
4070 “4175 
4023 4168 
4033 4170 
4045 1-4170 
-4048 1-4171 
4082 1-4172 


Rotating band columns 730 
Separating power data measured with 770 

770 
370 
370 
, : 180 
ures are presented in Fig. 9 as a func- 9.1 150 


the 6mm retating band column at 
constant rotational speed of about 


i rev. per sec under varying press- 


tion of boil-up rate. At atmospheric 9. &5 
pressures the separating efficiency , 1700 
proves to be fairly high (H.E.T.P. : 2900 
values of 15-25 mm being realized). 2150 
A tenfold reduction in pressure how- 2150 
1200 
1400 
640 


ever has a detrimental effect: a de- 
crease of about 70% in separating 
power is observed. 

Before discussing these pheno- 


_ 
— 


4123 1-4174 
“3921 1-4141 
“3935 1-4166 
“3951 1-4167 
3933 1-4167 
4008 1-4168 


260 
; : 2900 
attention should be paid to the rela- 2450 


mena more fundamentally, some 


tive vapour and liquid phase resist- ‘ , 2600 
ances under the conditions of the 3: 3 1500 


Vt Oo te & 


experiments. Considering the Rey- 8 370 


- & ts to te te « 


— ts oF me 


NOLDS numbers of vapour flow in 200 4048 “4170 


Tables 10 and 11 it should be con- f 2150 -3960 “4172 
2000 
2600 
390 
1000 
430 
470 
the n-heptane-methyleyclohexane sy- 170 


“3947 “4180 


] 

l 

cluded that the flow type is laminar. l 
“3945 1-4181 

1 

l 

l 

l 


Neglecting the influence of the rotat- 
4182 


“4184 


-4009 
“3975 
“3965 “4187 
4017 4192 
4108 1-4193 
-4142 1-4194 
about 2-7% of Rey, calculation from ‘ ‘ 230 1-3935 1-4154 
eq. (9) and (10) shows that (H.T.U.),; 36 150 )—-1-3952 1-4158 
presents only about 8°% of the total 36 9 14030 14152 
resistance to material transfer (assum- 36 ; 501-4088 1-4175 

230 1-3893 1-4135 
230—=«: 230 ~—-1-3900 1-4148 


ing element for a moment, the ma- 
terial transfer ii the vapour phase 
should obey eq. (9). At total reflux, 


tS to to ow ts ts te 
7 » _ 2 » : 
= 


l 
1 
1 
l 
I 
l 
l 
I 
l 
l 
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l 
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170 «6-1-4133 sd 4178 
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I 
I 
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I 
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stem liquid REYNOLDs numbers being 15 


ing liquid layer thickness at about 
0-1 mm [31]). From this discussion it 
may be concluded that at zero speed 
of rotation the overall efficiency may be adequately ing. However, the construction of the rotating band 
described in terms of vapour phase material transfer. column is such that the spinning band agitates the 
Now, at first sight, it might be expected that by in- liquid phase as well. In this study therefore the per- 
creasing the vapour phase material transfer by agita- centage liquid phase resistance will be assumed not to 
tion, liquid phase resistance would become controll- be influenced very much by agitation. Consequently, 
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Table 11. Performance of 6 mm rotating band column at reduced pressures 





Pressure Boil-up rate Rotational speed nD products Separating power 


Test mixture (H.T.U , 
—_ T.U Jor 
kN/m*? mm Hg mlh Rey rev sec Rep top bottom = 


n-decane-transdecalin : 1-4515 4672 
1-4376 -4672 

1-4327 4672 

1-4457 *4659 

1-4383 *4659 

1-4437 4659 

1-4396 -4650 

n-hexadecane- f 1-4584 -4880 
n-heptylbenzoate f 1-4663 -4906 
1-4644 4902 

1-4483 4903 





in the discussions to follow, vapour phase resistance of this paper, the data obtained should be plotted as 
to transfer will be considered controlling throughout. (H.T.U.)py/d Rey vs the ReyNoLps number of rota- 

In order to gain more insight into the mechanism tion Re,. The picture obtained by this method is 
of agitation in connection with the separation process, shown in Fig. 10. 


(ATU), 
d. Rev 
4 


THEORETICAL FOR 
ZERO ROTATION 
+ eg (9) 


Rey abt 220 


Rey abt 400 
1Rey abt 600 








4 6 8 1000 2 
REYNOLDS "NUMBER OF ROTATION 


WrdPy 
Lv 


Fig. 10. Influence of rotational speed at atmospheric pressure on overall vapour resistance in 6 mm spinning band column. 


a number of experiments were run at atmospheric It is seen that at Rep, up to 100 separating power 
pressure with n-heptane-methylcyclohexane under of the column equals approximately the theoretical 
wide variation of rotational speed in the 6mm column value calculated with eq. (9). Obviously the effect of 
(Table 10). As was suggested in the theoretical part rotation is negligible at the rotational speed applied 
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(about 2 to 3 rev. per sec). Increasing Re, above 100, 
a distinct influence of rotation is noted, however; it 
should be concluded therefore that Rep, about 100 re- 
presents the critical value discussed in the theoretical 
part. 


In the region Rep, 100-700 the data obtained at 
various vapour rates are well correlated by one single 
straight line. In this region the flowrate-separating 
power function is not affected by the stirring rate, 
(H.T.U.)joy/d Rey being constant at a specified rota- 
tional speed in agreement with the suggestion given 
in the theoretical part. 

When the Reynoips number of rotation is raised 
above 700, it has a less pronounced effect on causing 
In fact, 
220, rotation above Re, 


a further decrease in the (H.T.U.)oyp value. 
fora vapour flow of Re V 
1000 does not decrease (H.T.U.)oy any further. At 
higher vapour flow rates this critical Rep-number has 
a higher value. 


In discussing this phenomenon, the possibility of 
a limit set by liquid resistance which has so far been 
neglected should be considered here. From Fig. 10 
it is learned that at Rep = 1000, the overall (H.T.U.)oy 
has been reduced by about 75%. In the discussion on 
the relative values of phase resistances it was shown 
that total vapour resistance equalled about 92° when 
no agitation is applied; therefore when at Rep, > 1000 
liquid phase resistance should become controlling. 
a larger decrease than about 75% in (H.T.U.)o» would 
be expected ; moreover, no consideration has yet been 
given to the fact that liquid phase resistance is reduced 
too by the agitation. It seems therefore quite possible 
that liquid phase resistance as such cannot be the cause 
of the existence of a second critical rotational speed 
above which no further increase in overall separating 
power can be effected. 

When imagining the action of the rotating band 
on the separation process, it will be clear that the 
separating efficiency is particularly enhanced when 
the agitation increases the velocity of diffusion in the 
direction perpendicular to the directions of vapour 
and liquid flow. This effect will in particular be obtain- 
ed by the centrifugal action of the rotating band. 
However, agitation will increase diffusion in the direc- 
tion of vapour and liquid flow as well, causing a coun- 
ter-diffusional effect. Normally counter-diffusion or 
back-diffusion plays no important role [23], [32] in 
a wetted-wall column; at high rotating speeds how- 
ever the rotating strip becomes twisted and vertical 
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mixing effects in the liquid and vapour phases are 


becoming, very probable. 


Generally, back-diffusional effects are counter- 
balanced by increasing the convection streams in the 
separation process. Thus it should be expected that 
with higher liquid and vapour loads vertical mixing 
causes a smaller decrease in separating power than 
at lower loads. As already pointed out above, this 
effect is borne out by the experimental results: at 
lower loads (Rey = 220) separating efficiency is in- 
600, Fig. 10). 
The criterion whether the rotating band causes 


creased less than at high loads (Re, 


marked vertical mixing in vapour or liquid streams 
may be found in the relative linear velocities in the 
phase concerned and the circumferential velocity of 
the stirring mechanism. The higher the linear velocity 
in the phase, the less the probability that large por- 
tions of it will be swung back by the agitating element. 
Due to large difference in linear velocities, mixing 
effects in the vapour phase are not likely to be the 
same at atmospheric and at low pressure. Under both 
these conditions however in the experiments made 
liquid velocities are of the same order of magnitude. 
Therefore, when back-diffusional effects can be shown 
to be the same in the atmospheric and in the low- 
pressure experiments, liquid mixing should be regarded 
responsible. 

An indication of the back-diffusional effects is 
obtained by considering the (H.T.U.)oy-Rey relation- 
ship. When no vertical mixing occurs, (H.T.U.)oy is 
proportional to Rey. This is the case at low numbers 
of revolution, as shown in Fig. 11 (4 and 8-5 rev. per 
sec). At about 50 rev. per sec at atmospheric pressure, 
mixing causes (H.T.U.)oy to be proportional to Ref,”’. 
As shown in Fig. 11 too, this relation also holds at 
lower pressures. Now, recalling the fact that linear 
velocities in the vapour phase are 10-50 fold the 
atmospheric vapour velocities, it should be concluded 
that vertical mixing effects play only an important 
role in the liquid phase, thus limiting the ultimate 
separating power to be obtained by agitation. 

In connection with the back-diffusional effects 
found, it will be clear that the experiments at lower 
pressures can only be correlated with the atmospheric 
pressure experiments under conditions of correspond- 
ing vertical mixing, i.e. only experiments at the same 
rotational speed (50 rev. per sec) may be used in the 
correlation. 

The final correlation of the experiments with the 


6mm column at 50 rev. per sec is presented in Fig. 12 
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as curve A, in which the (H.T.U.)oy data, corrected 
for REYNOLDS number and ScuMipT numbers involved, 


(HTU)oy 
qd 


‘sorps | ; 
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ro LEGEND 
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Fig.11. Relationship between (H.T.U.)gyp and Reynoips 
number of vapour flow at 
(6mm spinning band column) 


various rotational speeds. 
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are plotted against the ReEyNoLDs number of rotation. 
It is seen that the data obtained at various pressures 
without varying the rotational speed correlate in a 
like manner as data obtained at atmospheric pressures 
with variation of rotational speed (see curve B). Ob- 
viously correlation of data against the ReyNo.ps 
number of rotation is sound; it shows that an increase 
in density has the same effect as an increase in rota- 
tional speed. At pressures of about 2-0 kN/m? (15 mm 
Hg) the density of the vapour phase is so low that even 
agitation at 50 rev. per sec does not have any effect 
on the overal separating power. As may be calculated 
from eq. (9), separating power is even 30% too low, 
owing to vertical mixing effects in the liquid phase 
at this high number of revolution. 


In Fig. 12 the results obtained with the 3 mm diam- 
eter spinning band column are plotted too (curve C). 
Further data obtained at atmospheric pressure by 
Brrcu [3] et al. on a 36mm diameter rotating band 
column are included (curve D). In both cases an in- 
fluence of REYNoLDs number of rotation on the sepa- 
rating power analogous to that found with the 6mm 
diameter column is shown. At rotational speeds of 
the Brrcn column higher than that corresponding 


50 REV P SEC 
VERTICAL MIXING IN 
LIQUID PHASE 


_-THEORETICAL FOR ZERO AGITATION 
~~ eq (9) 


2-20 REV. P SEC. 
NO VERTICAL MIXING 
IN LIQUID PHASE 





7000 70000 
Rep. REYNOLD S'NUMBER OF ROTATION wr dpy 
By 


Fig. 12. Correlation separating power data of spinning band columns 
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with Re, ~ 6000, (H.T.U.joy is found to increase 
again, showing the same mixing effects as described 
with the 6mm column. 

In spite of correlation against Re, it is seen from 
Fig. 12 that the data for the three columns do not 
coincide on a single line, but three separate functions 
are obtained. This should be ascribed to different 
(Rep) 
might be expected on the ground of the theoretical 


ert, Values for each of the three columns, as 


discussions on agitation. As suggested, the critical 
values for Rep, should be a funetion of the ratio of 
diameter to clearance of the columns concerned. The 
following (Rep)... 


tabulated from Fig. 12 and the dimensional data on 


values and dA values may be 


the three columns: 





Diameter (Re p).., 


adiA 
(appror.) 


column 


(Re p) 


mm 


1-7 
1-7 





It is seen that (Rep) increases markedly with 
the dA ratio. 


Actually, it may be calculated that Re, is pro- 


crit 


portional to d/A raised to the 2-5 power. The influence 
of column dimensions in this respect is thus much 
more pronounced with spinning band columns than 
with concentric tube columns [e/. eq. (11)]. This quan- 
titative relationship found should be considered with 
some reserve however; it is only indicative of the 
trend which various factors involved have on the 
efficiency of the rotating band column. 

Under conditions liquid when mixing effects play 
no role, the following ultimate relation for (H.T.U.oy 
may be derived from Fig. 12. 


(H.T.USoer ~ 
d 


0-065 | : ) "Re nr. Rey . Sc}". 


Though this equation has only an illustrative value, 
it shows that at constant boil-up rate, Rey and Sey 
being approximately independent of total pressure, 
separating power is controlled by Rep. Thus, at con- 
stant speed of rotation, vapour density being pro- 
portional to pressure, separating power varies with 
pressure, raised to approximately 0-6 power. 

The results gained with the spinning band columns 
may be compared with those obtained in other transfer 
processes influenced by agitation. 
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WILLINGHAM [35] and co-workers gave data on a 
rotating concentric tube column. Their results clearly 
show the existence of a critical rotational speed, above 
which a transition range follows in which the separat- 
ing power increases sharply with rotational speed; at 
higher numbers of revolutions separating power in- 
creases less markedly. Though the results do not cor- 
relate well when using the method of the present paper, 
in the higher ranges of rotational speed (H.T.U.)oy 
seems to decrease Re',”. 

Measurements on material and heat transfer in 
agitated vessels were made by Hixson [19], [20] et al. 
and by Curtton, Drew and Jesens [10]. Using padd- 
le type stirrers, they found the material and heat trans- 
fer coefficients to increase with Rep raised to 0-62-0-67 
power. 

In both cases of the rotary concentric tube column 
and the agitated vessels material and heat transfer 
are enhanced to about the same degree as found with 
the spinning band columns. It therefore seems justified 
to consider the spinning band column as a series of 
agitated units, in each of which the effect of agitation 
has the same result as found in separate agitated 
vessels. 

Finally it should be concluded from the present 
study that vapour-phase agitated columns cannot, as 
a rule, be regarded as special low-pressure columns. 
In fact, at pressures of about 1-3 kN/m? (10 mm Hg) 
and below, the separating efficiency is not better than 
that found with a wetted-wall column of the same 
dimensions and length. 
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NOTATION 


a = constant in eq. (1) and (2) 
constant in eq. (11) 
diffusivity 
- diameter of column m 
height of transfer unit m 
liquid flow rate kmols/sec 


m = averaged slope of equilibrium line 
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(see footnote on page 175) 

theoretical plate number of column 

constant in eq. (1) and (2) 

constant in eq. (1) and (2) 

ReyNoips number of liquid flow, 4wuyo,/“, 
ReyNoips number of rotation, orde,/uy 
REYNOLDS number of vapour flow, dupoyp/ Hy 
half width of rotating strip m 
Scumipt number, u«/Do 

temperature K 
average linear velocity m/sec 
vapour flow rate kmols/ sec 
thickness liquid layer m 
liquid composition, mol fraction 

liquid composition, in equilibrium with vapour 
vapour composition, mol fraction 

vapour composition, in equilibrium with liquid 
height of column m 
clearance bet ween two rotating cylinders or between 
column diameter and spinning band, (d—2r) m 
viscosity Nsec/m? 


density kg/m?* 


angular velocity rad/sec 


Indice 


lower end of column 

upper end of column 

liquid phase 

vapour phase 

at plane of contact between liquid and vapour 


phase 
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The flow of compressible fluids in pipes is governed by 
two relations which can be expressed mathematically 
as the energy balance or Bernoulli equation and the 
equation of state of the fluid. The energy balance is 
represented by the equation 

V?—V? 
Z,) + (P22 — Pi %) + (| 5-*] + 


29e 
+ J (ug —u)=JQO+N, 


(Z. 


or in the differential form 


, VaV , 
dZ +d(pv) + * Jdu=JdQ+dW,. (2) 

When a fluid is flowing in a pipe or nozzle the ex- 
ternal work term W, is zero, but friction causes a de- 
gradation of mechanical energy into heat and from 


the first law of thermodynamics 
JdQ+dF =Jdu-+ pdv. (3) 


The frictional loss in a pipe is usually expressed by the 
equation 
4/V'dL 


dF 29. D 


(4) 
and from eqs. (2), (3) and (4) the mechanical energy 
balance is obtained: 

Vav | 4fV dL 


dZ +-vdp-+ ge : 29D 


0. (5) 

The friction factor / is a function of the Reynolds 
number and of the roughness of the pipe surface. The 
Reynolds number for the isothermal flow of a perfect 
gas in a uniform pipe is constant along the pipe since 
the mass flow rate G = o V is constant and the viscosi- 
ty of a perfect gas is independent of its pressure. Any 
change in the temperature of the gas will alter its visco- 
sity and therefore the Reynolds number and friction 
factor, but for fully turbulent flow (Reynolds numbers 
greater than 10000) the change in the friction factor 
is relatively small and may usually be ignored. 

Eq. (5) may be written, substituting V — Gv, as 


dZ +vdp+ —o0e + oom 


= + o,p eel =o. (6) 
c c 


This equation can be integrated only if the relation 
between p and v is known. For isothermal flow pv 
constant, and since the term dZ can usually be neglect- 
ed when dealing with gases, eq. (6) may be rearranged 
” @ dv . 4/@ 

P 2) vv, 7 - 
4 . } L=0. 
Pay ge e 29-D oe @) 

Integrating between state | and state 2 for a pipe 
of length L gives 


L 
ef 
- p|dl=o 


2 
l @ fdv 
Pai’i. 
1 


; 
| pdp + 


pv, @ 


|2In 
Ge(Py + Pe) 


Py 





29e Gav 
where o,, is the average density of the gas at the 
average pressure (p,; + p,)/2. When the term 4/L/D is 
large in comparison with 2 /n(v,/r,), so that the latter 
may be ignored, eq. (9) is identical to the friction fac- 
tor equation for an incompressible fluid, with the gas 
density taken as the average between the initial and 
final states. 

The equation of state for the adiabatic flow of a 
perfect gas, derived from the thermodynamic rela- 
tions, is 

x—1 VdV 


d(pv) + ~ = 


= (10 


= constant 
= Py 7 
dv 
vdp=d(pv) —(pv) 
1 (* -1\ @ 
2 me Ge 


x—1 ) P| de 
29c { v 


vdv 





= {Pat + 
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Substituting in eq. (6) and dividing by v* gives 


| x ] 


. Gtvt) dv 
) Pa ") 


294 { 73 


| x—I 


as : tts 


2g f \r3 


13 
jl 4/L a | M9) 
ley aj" D : 


-I\ oe 


| 
x 2 | 


The flow conditions must satisfy the equation of 
state, eq.(11), and the energy balance, eq. (13), but 
owing to the complicated nature of these equations, 
a direct mathematical solution is out of the question 
and a numerical trial-and-error method must be foll- 
owed, For example, if it was required to calculate the 
pressure drop for a known flow rate @ in a pipe of 
specified dimensions, it would first be necessary to 
assume a value for the specific volume of the gas at the 
exit of the pipe, v,, and substituting the known values 
of py, 1%, G, D and f, to calculate the pressure at the 
exit, py, from eq.(11), and then the length of pipe L, 
necessary to give this pressure drop using eq. (13). 


Then if the caleulated pipe length differed from that 


values of v,, and repeat the calculation until the cor- 
rect conditions were found. 

A convenient graphical solution to these equations 
has been developed in terms of three dimensionless 
ratios, (py — Pe) py, 4/L,D and the initial Mach number 
Ma,, where 


Ma, - 
Voy 


Gv, (14) 
| xg Mat 
Substituting in eq. (13) 
1; 1 
x | Maj 


Also from eq. (11) 


Pils P22 (v3 


296 
(x —1) Ma? % \7)) . 
et (t— (28) |}. (16) 


The chart shows the pressure drop ratio (py — ps/Py 
as a function of the initial Mach number Ma, and the 


friction number 4/L/D for adiabatic flow when the 
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ratio of the specific heats x is 1-4. The pressure drop 


ratio has also been calculated for isothermal flow 
using eq. (9) which may be written as 
@ 


29e Pi Oar 


D | 
4/L) 
D {° 


. 
2In 34 


an Ma; —; 


» 
- Var 


12 In 
| 
For values of Ma, below 0-25 the calculated pressure 
drops for adiabatic and isothermal flow are almost 
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Fig. 1. 


identical and are represented by single lines on the 
chart. Appreciable deviations occur between adiabatic 
and isothermal flow in the top right hand corner of the 
chart, when the full curves represent adiabatic flow 
and the broken curves the isothermal state. 

For each value of 4/1, D there is a critical pressure 
drop at which the velocity at the discharge end of 
the pipe reaches the velocity of sound in the gas. At 
this point the discharge rate and initial Mach number 
Ma, attain their limiting values. The critical pressure 
drop and the limiting value of Ma, for adiabatic flow 
are indicated by the broken curve across the top of the 
chart. 

The chart may be used, in conjunction with 
eq. (14), to caleulate the pressure drop in a pipe for 
a given rate of flow, or the rate of flow when the pressu- 
re drop and dimensions of the pipe are known. The 
resistance due to pipe bends and fittings may be 
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evaluated in terms of the equivalent length of straight 


pipe, or in terms of the number of “velocity heads” 


lost, as recommended by Lappe [1]. 
The friction number 4/ LD represents the number 
of “velocity heads” lost through friction, since from 


the friction factor equation 


(o V*/29¢) 
where 9 V?/2g. represents the “velocity head” o1 
pressure equivalent to the kinetic energy of the fluid. 

When the chart is used to calculate the rate of 
flow of gas through a pipe connected to a vessel, the 
initial pressure p, must be taken as the pressure within 
the pipe at the inlet, this pressure being lower than the 
pressure p, in the vessel by an amount equivalent to 
the velocity head of the gas, plus any frictional losses 
at the inlet. In the case of a rounded inlet the frictional 
losses are small and the pressure drop is equivalent to 
the kinetic energy gained by the gas. Since this is 
a reversible process the pressure drop at the pipe 
inlet should be evaluated separately and added to 
the irreversible pressure loss due to friction in the pipe, 
to obtain the total drop from the receiver to the pipe 
outlet. 

However, for many purposes, a sufficiently accura- 
te result is obtained by assuming that the pressure 
drop at the inlet is an irreversible loss, and adding 1-0 
velocity head for a rounded inlet, or 1-5 for an abrupt 
inlet, to the friction number 4/ L/D. 

An exact method for calculating the pressure drop 
for a gas flowing from a vessel through a pipe with 
a rounded inlet has been presented by Larrie [1]. 
A detailed comparison of results calculated by Lapp- 
LEs method and using the present chart with an allow- 
ance of 1-0 velocity head for the inlet shows that the 
deviation is appreciable only when the friction number 
for the pipe is less than 2-0. When the friction number 
for the pipe is 1-0 the maximum discharge rate calcul- 
ated by LarrLes method is 4% higher than that calcul- 
ated from this chart, assuming a loss of | velocity head 


at the inlet. 


Example 
Calculate the rate of discharge of air from a reservoir 
at 150 lb in? gauge and 70° F through 33 ft of pipe, 
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2-067 in. internal diameter, with an abrupt inlet and 


three elbows. 


Assumed friction factor, f = 0-0040 
Friction number, 4/L/d 
Pipe = 3-06 
Inlet (equivalent) = 1-5 
Elbows (equivalent) = 1-44 


Total 6-0 
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0-0233 ft?. 


6-49 lb’sec. 


For the same problem the rate of flow calculated by 
LapPcLr [1] (page 423), was 6-41 lb see. 


List of symbols 
Pipe diameter 
Friction factor 
Gravitational constant 
Mass flow rate 
Mechanical equivalent of heat 
Pipe length 
Mach number at inlet of pipe 
Pressure 
Heat added to system 
Work done on system 
Internal Energy 
Specific volume 
Velocity 
Velocity of sound 
Height above a datum plane 
Gas density 


Ratio of specific heats (¢,/¢,) 
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